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Abstract 

This investigation is concerned with the dynamics of a gravity- 
gradient stabilized flexible satellite in the neighborhood of a deformed 
equilibrium configuration. First the equilibrium configuration is deter- 
mined by solving a set of nonlinear differential equations. Then stability 
of motion about the deformed equilibrium is tested by means of the 
Liapunov direct method and the natural frequencies of oscillation of the 
complete structure calculated. The analysis is applicable to the RAE/B 
satellite. 
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1. Introduction 


With the advent of large spacecraft, flexibility has become an in- 
creasingly important factor in the system attitude stability. Early de- 
signs of spacecraft were based on rigid body analysis, according to which 
rotational motion is stable if it takes place about the axis of maximum or 
minimum moment of inertia and unstable if the body rotates about the axis 
of intermediate moment of inertia (see, for example, Ref. 1, Sec. 6.7). 

The erratic behavior of the Explorer I, a satellite stabilized about the 
axis of minimum moment of inertia, prompted a re-examination of the rigid 
body assumption. Indeed, Thomson and Reiter^ were able to attribute the 
behavior of the Explorer I satellite to energy dissipation resulting from 
the vibration of flexible antennas. This conclusion was corroborated by 
Meirovitch. 3 References 2 and 3 used the so-called "energy sink" approach. 
Their main conclusion was that a flexible satellite cannot be stabilized 
about the axis of minimum moment of inertia, leaving as stability criterion 
what has come to be known as the "greatest moment of inertia" requirement. 

For a number of years, no significant additional work on the 
stability of flexible spacecraft was performed. Some work on cable-connected 
space stations cannot be really considered pertinent. Some investigation 
that can be regarded as being related to flexible spacecraft is that by 
Hooker and Margulies,^ who model a satellite as "a set of n rigid bodies 
interconnected by dissipative elastic joints," and forming so-called 
"topological trees." 

The first serious attempt to treat rigorously the flexibility ef- 
fects on the attitude stability of flexible satellites can be attributed 
to Meirovitch and Nelson. 5 Reference 5 investigated a satellite with 
elastic appendages by means of an infinitesimal analysis. It appears that 
Ref. 5 uses modal analysis for the first time in conjunction with the 
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stability of flexible spacecraft. At the same time. Nelson and Meirovitch 
used the Liapunov direct method to investigate the stability of a rigid 
satellite with elastically connected moving parts. The motion of a satel- 
lite consisting of two rigid bodies connected by an elastic structure was 
investigated by Robe and Kane 7 by means of an infinitesimal analysis. Simu- 
lating a spacecraft by a set of rigid masses interconnected by massless 
elastic members, Likins 8 derived the corresponding equations of motion, and 
indicated that a solution can be obtained by modal analysis. Reference 8, 
however, does not produce an algorithm for the solution of the equations. 

9 

Thermal effects and solar radiation pressure were found by Etkin and Hughes 
to cause the anomalous behavior of spinning satellites with long flexible 
antennas. The flexibility effects on the attitude motion of spacecraft 
were also investigated by Modi and Berenton^ 8 but the validity of their 
analysis is in doubt, as they restrict the satellite vibration to planar. 

An interesting paper by Newton and Farrell^ presents a method for 
evaluating the natural frequencies of a flexible gravity-gradient stabi- 
lized satellite. In the process. Reference 11 linearizes the equations of 
motion about the deformed equilibrium. As generalized coordinates, the 
investigators consider complete deformation patterns of the satellite. 

This procedure is not only unorthodox but also tends to limit the number of 
degrees of freedom of the simulation, not to mention the fact that one must 
guess in advance configuration patterns. Moreover, there is some question 

as to the evaluation of the equilibrium configuration. Nevertheless, the 

12 

paper contains some interesting ideas. A paper by Likins and Wirschmg 
proposes to introduce the concept of "synthetic modes" in conjunction with 
a "hybrid" coordinate system, where the latter is defined as a set of co- 
ordinates consisting of rotational coordinates of the spacecraft as a whole 
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and modal coordinates for the flexible appendages. This idea, however, 
was introduced earlier in Reference 5. 

All preceding investigations have one thing in common, namely, 
they are all discretization schemes. Some use lumping of the distributed 
mass of the elastic members, a procedure referred to sometimes as spatial 
discretization, and others use series truncation in conjunction with modal 
analysis. In a first attempt to apply Liapunov's direct method to hybrid 
systems from the area of satellite dynamics, i.e., without using any dis- 
cretization scheme, Meirovitch 1 ^ 1 ^* 1 ^ studied the stability of spinning 
rigid bodies with elastic appendages. It should be pointed out that the 
term "hybrid" refers here to a system defined by sets of both ordinary and 
partial differential equations, a concept different from that used by 
Likins and Wirsching. 12 Several new ideas were introduced in Ref. 13, such 
as the use of the bounding properties of Rayleigh's quotient to eliminate 
spatial derivatives from the problem formulation and the use of testing 
density functions. 

The ideas of Refs. 13-15 have been pursued by Meirovitch and 
Calico 16 ’ 17 for the case in which testing density functions cannot be de- 
fined readily. References 16 and 17 develop the so-called "method of 
integral coordinates," whereby certain integrals are identified and defined 
as generalized coordinates. Then, using the bounding properties of 
Rayleigh's quotient as well as certain Schwarz's inequalities for functions, 
a function x bounding the Hamiltonian H from below is obtained, k <_ H, so 
that x can be used as a testing function in conjunction with Liapunov's 
direct method. The method of integral coordinates is basically a dis- 
cretization scheme. 

One problem that has received little attention in the technical 
literature is that of deformed equilibrium, which can be referred to 
mathematically as "nontrivial equilibrium." Such problems arise in the 
case of gravity-gradient or spin-stabilized satellites with very flexible 
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appendages that are not aligned with the satellite's principal axes. Find- 
ing the equilibrium configuration can be quite a problem in itself, 
especially if the governing equations are nonlinear. Addressing himself 
to this problem, Flatley 18 obtained the nonlinear equilibrium configuration 
of the Radio Astronomy Explorer (RAE) satellite by means of an analogue 
computer. Deformed equilibrium has also been considered in Ref. 11, but 

the details are not clear and no plot of the deformed equilibrium is shown. 

19 

In seeking stability statements for the RAE/B satellite, Meirovitch ob- 
tained as a by-product the nonlinear deformed equilibrium, thus confirming 
the results of Ref. 18. 

The present study is concerned with the stability of a hybrid dy- 
namical system about nontrivial equilibrium. It contains many of the for- 
mulations and results of Ref. 19. Qualitative stability statements are ob- 
tained for the RAE/B satellite by both the Liapunov direct method and by an 
infinitesimal analysis. In connection with the infinitesimal analysis, the 
natural frequencies of oscillation about the nonlinear nontrivial equilib- 
rium were obtained by a method developed by the first author of this re- 
port. 20 The method of Ref. 20 considers a state vector consisting of 
generalized coordinates and velocities, where the coordinates include both 
rotations and elastic deformations , and develops an eigenvalue problem in 
terms of real quantities alone. The stability statements of Ref. 19 and 
corresponding statements obtained from the solution of the eigenvalue prob- 
lem agree completely. 

2. Problem Formulation 

We shall be concerned with the motion of a body consisting of n + 1 
parts, of which one part is rigid and n parts are elastic. The domain of 
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extension of the rigid part is denoted by D 0 and those of the elastic parts 
when in undeformed state by (i = l,2,...,n) (see Fig. 1). Correspond- 
ingly, the masses associated with the domains D-j are denoted by m i 
(i = 0,1,..., n), so that the total mass is m = . 1 ^ m^. The elastic domains 

are rigidly attached to Dg and have common boundaries only with Dq. 

The body m is assumed to move in a central -force gravitational 
field, with its mass center describing a given orbit about the center of 
force C.F., where the latter is assumed to be fixed in an inertial space. 

In describing the motion of m it will prove convenient to identify 
a system of axes xyz (see Fig. 1) with the undeformed state. The origin c 
of xyz is taken to coincide with the mass center of m in the undeformed 
state and axes xyz themselves coincide with the principal axes of m in the 
same state. Note that the system xyz is embedded in the rigid part Dq but 
is not necessarily a set of principal axes for that part. We shall assume 
here that the nature of the elastic motion is such that the mass center of 
the entire system remains at the origin of xyz. In measuring elastic de- 
formations, we consider reference frames x^zj fixed relative to the 
elastic domains D i (i = l,2,...,n), where the direction of these axes is 
chosen parallel to that of the elastic deformations. The origin of axes 
x-v-z- is denoted by 0.: and in general it need not coincide with c. 

I* 7 ! 1 1 

Next let us denote the radius vector from the mass center c to a 
point in the domain (i = 0,1, ...,n) by hj + r-j , where the point coincides 
with the position of an element of mass dm^ when the body is in undeformed 


state. The constant-magnitude vector h;j denotes the radius vector from c 

to 0-j ; clearly h 0 = Q.- 0n the other hand > li 1s the radius vector from °i 

to the point in question, and its components represent the independent 

spatial variables associated with a point in the domain D-j . Denoting by ij , 

j. and k-; the unit vectors along axes x.j , y-j and z^ , respectively, we can 
' ■ V '' 

write h.j + r n - = (h x1 * + x.j)iLj + (hyj + y n *) j-j + (h 21 - + )k n - (i - 0,1,..., n). 
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In describing the elastic deformations, we can use the Lagrangian or the 
Eulerian approach. According to the Lagrangian approach the independent 
variables are those of the body in undeformed state, whereas in the 
Eulerian description of motion the independent variables are those of the 
body in deformed shape. For infinitesimally small deformations the two 
points of view coalesce, but for large deformations they do not. When it 
is necessary to calculate the stresses in a body undergoing large deforma- 
tions, the Eulerian approach is more convenient. Although we shall be con- 
cerned with relatively large deformations, we have no interest in the in- 
ternal stress distribution, and because of kinematical considerations we 
shall find it more convenient to use the Lagrangian approach. Hence, de- 


noting by u-j the elastic displacement vector of dm^ , and recognizing that 
the vector depends both on spatial position and time, we can write it in 
the form £j = u., (x^ ,y^ ,z n - ,t)i^ + V n - (x n * ,y n ‘ ,z^ ,t) + w^ (x-j ,y^ ,z^ »t)k-j , where 
u-j , v-j and w i are displacement components measured along x-j, yj and z n - , 
respectively. If R c is the radius vector from the center of force C.F. to 
the mass center c, then the position relative to the inertial space of a 
mass element dm i in deformed state is given by R^-j - R c + h-j + + Uj * 

It should be noted that, by the definition of the mass center, 

i=0 / mi (ti + li + Hi> dm 1 = .9- 

In view of the above discussion, the kinetic energy can be written as 


T = H 5di • B 


n 


f 


di dm i = l m Sc • Be + i ^ k k + ii 


+ Mi) • ( h li + ii + Hi ) dm i 


( 1 ) 


where the first term on the right side of Eq. (1) is recognized as the kinetic 
energy of translation of the mass center c and the second one as the kinetic 


6 



energy due to motion relative to c. Dots. denote derivatives with respect to 
time. Denoting by w the angular velocity of the set of axes xyz, hence also 
of the sets x^z. (i = l,2,...,n), and recalling the expression for the 
time derivative of a vector expressed in terms of rotating coordinates, we 
obtain 


fii + E, + Hi = Hi + & + IT + Hi > 

in which u! = i - + vj-r + wk- is the velocity of dm i relative to c due 

M | I iU | 1^1 ^ I 

to elastic effects alone. Introducing Eq. (2) into (1), we arrive at 

T = r ■ ic • 5c + 1 a • id • » + - • I J ra . % + n + hi) x Hi dm i 


( 2 ) 


\ f Hi • Hi dl "i 

1=1 •’m,* 


(3) 


where Jj is the inertia dyadic of the body in deformed state about axes 
xyz. 

Equation (3) is most conveniently expressed in matrix form. The 

matrix forms of the vectors R„ and u are simply (R c ) and M , respectively. 

The inertia dyadic and the term on the right side of Eq. (3) require 

~ . n (i) 

further attention. The inertia dyadic can be written as > 

where J W (i = 0,1, ...,n) is the inertia dyadic associated with domain 

D. when the corresponding mass is in deformed shape. The superscript 

i indicates that the dyadic is expressed in terms of the base x^y^z^. 

This would require that we express w in the same base. It is simpler, 

however, to express every in the base xyz instead. Denotinq the 

vector r- by rf^ and r! 1 ^ when expressed in the base xyz and x.y z , 

V-1 «*1 ~] . Ill 

respectively, and by {r? 0 ^ > and (r. } the associated column matrices. 
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the relation between the two can be written as {rj } 3 [£-,*] {rj" }, 
where [#,•] is the matrix of direction cosines between axes x^y ^ and 
xyz. In a similar fashion, if we denote by 0^ and the inertia 
dyadics when expressed in the base xyz and x.y.z. , respectively, and 


by [J^] and [J^V] the associated inertia matrices, then the relation 
between the two can be shown to have the form = C‘i 3 T C J dl 3CM3- 

With the understanding that the inertia matrices imply the body in deformed 
shape, we can drop the subscript d. Moreover, we shall drop the super- 
script i when it agrees with the subscript. Hence, the inertia matrix for 
the entire body, expressed in the base xyz, takes the form [J^] = 

.| 0 [^•] T [J i ](> i ]. We note that O 0 ] = [1], where [1] is the unit 
matrix. A similar analysis can be performed with regard to the third 
term on the right side of Eq. (3). It follows that Eq. (3) can be written 
in the matrix form 

T • 1 " T _ _T 


i(1) 


T = ^ m (R c 


t - < [*<"> 


* * 
+ r. + 


(4) 


} {R c > + ~ [>.■] [0 ^ ][«- i + M i I 

i=0 i=l J m 1 

(0) T- 1 n f • T • 

u. ][*,-] {uijdm., + 5- z {u : } (u'-}dm. 

1 1 1 c i=1 Jm i 1 1 

where [h^ ^ + u- b is a skew- symmetric matrix whose elements satisfy 

the relation h<°> + rJJ + uffi = j, < h U + r U + 1" which 


nnw 


is the epsilon symbol (see Ref. 1, p. 109). Clearly, {ul> represents 


the matrix notation of u ! . < 

The potential energy results from two sources, namely, gravity and 
elastic deformations, denoted by Vq and Vg, respectively, so that 
y = Vq + V^. From Ref. 15, we conclude that the gravitational potential 
energy can be written as 
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where tr denotes the trace of a matrix, and U a > is the column matrix of 
direction cosines between the direction of the vector R c and axes xyz. 

The elastic potential energy, also known as strain energy, requires 
special attention, particularly in the case of large deformations. No 
general expression, such as for T and Vq, can be written for Vp. This is 
so because an explicit form requires the knowledge of the type of elastic 
members involved. For the moment, we shall be content to write 

n 

V c = Z V r . (6) 

11 i=l 

where (i * l,2,...,ri) is the elastic potential energy associated with 
the member occupying the domain when the member is undeformed. We 
shall return to the elastic potential energy shortly. 

At this point it appears desirable to determine the functional 
dependence of the kinetic and potential energy in order to derive general 
Lagrange's equations of motion. To this end, we must specify the nature of 
the elastic members. We shall be concerned with one-dimensional members 
capable of flexure in two orthogonal directions. Any axial displacements 
will be assumed to be a result of change of length caused by the trans- 
verse displacements and not because of axial flexibility. In essence, the 
members are cantilevered bars undergoing large transverse displacements 
(see Fig. 2). Although we shall use nonlinear theory for the elastic 
motion, this will be because geometric nonlinearities and not as a result 
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of nonlinear stress-strain relations. The mass distribution is arbitrary. 


but some of the members carry tip masses. 

Letting the radius vector r^ be aligned with axis when the bar 
is undeformed, we conclude from Fig. 2 that 



x i 


i . 


i = l,2,...,n 


and 



■ “i ( x i ’t) = + w i (x i ,t)k i 



( 7 ) 

( 8 ) 


In view of this, the elements of the inertia matrix for the rigid member 
can be written as 


°on ■ A 0 ’ J 022 = B 0 ’ J 033 = c 0 


J 012 = °021 ' °013 = J 031 * J 023 ‘ J 032 0 


(9) 


where Aq, Bq, Cq are the principal moments of inertia of the rigid part, 
whereas these for member i are 


°m 


J i22 


J i33 


J il2 


J il3 


1 p i [(h yi + v n -) 2 + (h zi + w i ) 2 ]dx i + m i [(h yi + v i ) 2 + (h zi + w i ) 2 ] 


f l 1 

0 


P i [( h xi + x,) 2 + (h zi + w i ) 2 ]dx. + nii C(h xi + x i ) 2 + (h zi + w.) 2 ] 


f \[(h xi + x i ) Z + (h yl - + v i ) 2 ]dx. + + x,) 2 + {h yi + v,-) 2 ] 

J i21 = ‘ L 1p i (h xi + x i )(h yi + v T )dx i ' m i (h xi + x i )(h yi + v 1> 


J i31 


■Hi 

-0 


Pl ( h xi + xpth^. + w.)dx. - ffli (h xi + x i )(h z1 + wp 


X i = *i 


x r £ i 




Note that and m. are mass densities and tip masses, respectively, and 

h . , h . , h . denote the coordinates of the points of attachment of the 
xi yi zi 

booms measured from the mass center along axes x^y^z^ (i = l,2»...,n). 

We shall not specify the mass densities at this point. 

The desired equilibrium configuration is that of gravity-gradient 
stabilization. That implies that the mass center c moves in a circular 
orbit with the constant angular velocity (see Fig. 3), and the set of 
axes xyz coincides with a set of orbital axes abc, where a^ coincides 
with the direction of the radius vector R , b is tangent to the orbit 
and in the direction of the motion, and c is normal to the motion. Note 
that the orbital axes abc rotate relative to an inertial space with angular 
velocity Jl about axis c. The orientation of axes xyz with respect to abc 
is given by three angles 6j and {w} depends on these angles and angular 
velocities 6j (j = 1,2,3). Because the first term in the kinetic energy, 

Eq. (4), is constant for a circular orbit, it will be ignored in future 
discussions. Moreover, the last term depends on the elastic velocities, 
so that the functional dependence of T is 

T = T(6j, e j , v . , v i , w i , w. ) , j. = 1 ,2,3 ; i = 1 ,2,. . . ,n (11) 

The gravitational potential energy V G contains the matrix {£ a >, which 
is defined as the matrix of direction cosines between R c and xyz. Since 
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xyz can be obtained from abc by means of the rotations 0j (j - 1*2,3), it 
follows that 


Vq - Vg (Qj, » w. ) , j = 1,2,3 *, i - l,2,...,n 


( 12 ) 


It remains to establish the functional dependence of V E . This 
requires some elaboration, particularly because of the geometric nonlineari- 
ties involved. First we wish to distinguish between the potential energy 
V EA due to axial motion, and the potential energy V EB .due to flexure. 

Next let us consider Fig. 4 and denote by s i the distance to any element 
of mass dm.j when measured along the deflected bar and by x^ when measured 
along the original direction of the undeflected bar. We shall assume that 
the bar is inextensional , so that these two distances remain the same, 
s . = x _j . An element of length along the deflected bar can be obtained 

from 

(ds^ 2 = {dx i + di^) 2 + (dv t ) 2 + (dw^ 2 (13) 

Assuming that du. is one order of magnitude smaller than dv^ and dw^ , 

recalling that ds i - dx - , and rearranging Eq. (13), we arrive at 

1 rr dv il 2 f dw i] 2 l . . I , o 

du i = - \ [jd^j + [d)qj J dx i ’ 1 = 


(14) 


so that the axial displacement resulting from the transverse displacements 
is negative. Because for inextensional motion the axial force P x ^ does 
not depend on the axial displacement, and, moreover, because a tensile 
force opposes the motion, we have 


V 


EA 


n 

i-l i 


£ i 


P xi du i = 2 * 


r£ i p 1 

rm 2 . i 

j ' 

t 2 i 

p xi 

J n 

.UxiJ ' 


1 J 


dx i 


(15) 
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where total derivatives have been replaced by partial derivatives in 
recognition of the fact that the displacements are functions not only of 

spatial position but also of time. 

The flexure potential energy is due to the bending displacements v^ 
and w r We shall denote the bending moment associated with the displace- 
ment v 1 by M zi and the change in slope corresponding to an element of 
length in deformed state by d* zi because they both take place about the 
z.j-axis. Accordingly, the analogous quantities associated with flexure 
about y i are denoted by M yi and d* yi , respectively. It follows that the 
flexure potential energy can be written as 



" f 


i=l 



06 ) 


But the bending moments M y ^ and M Z1 - can be written in terms of the 
associated flexural stiffness and radii of curvature, as follows 



El 


J£1 


yi 




where, from Fig. 5, the radii of curvature have the form 


(17) 


ds . ds . 

d = yi r = 21 

K yi " d<j> yi J Zi d* 2i 


(18) 


in which 


ds - 
zi 


1 + 


dv^2 


r qv ii 

U x i. 


1/2 


(dw^ 2-i 1/2 dx 

dx i . ds yi = L 1 + W J 1 


09 ) 


Moreover 


dv-j ^ 

5T = tan ♦zi 


dw . 

; ta " V 


( 20 ) 
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From Eqs . (20), it follows that 


-1 dv-; 

dxi 


d<j> zi - = d [tan 


fd^' 

< 


TT 

dx- 


dx. 


% 5 d ( tan dx. 


-1 dw-j'j 
i 


1 + 


dv,*i 2 




fdw^ 




l dx ii 


dx^ 


1 + 


dw.p 2 


dx i-» 


1 + f. 


d 2 w i 

dx-^ 

7 dx i 

dw,' z 1 


idx-j 


(21) 


Finally, introducing Eqs. (17) through (21) into Eq. (16), we obtain 

2 


I n f 

EB " 2 1*1 J M 


1 


n 


Vcd - „ F I (^zi d< ^zi + ^yi d ^yi^ ” 2 ^ 


El 


(d* z1 )‘ 


zi ds 


Zl 


(%•) 

+ EI yi 


2 i 


Li_ _I " f z i 

' j i=l Jo 


( d V 

1 2 

, d *i 2 

i 


r-T V> i J + FT 

21 r,.r dv il2l5/2 ^ 


r r u v 1 



^d^w. 

idx-j 2 . 

2 * 

. . 

[ 1+ l 

f dw r 

1 dx * 
L ax lJ 

,2 

i 1 

CJ1 

INJ 


}dxi 


(22) 


Recalling that the flexural displacements depend also on time, and writing 
binomial expansions for the denominators in Eq. (22), we arrive at 
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where the terms involving av-j/ax^ and aw^/ax 1 - are recognized as the cor- 
rections due to the geometric nonlinear effect. 

In view of the above, the potential energy has the general functional 


form 


V E ' V EA + V EB 


Vr (v! , v" , w! , w 1 .') 

1 1 1 i 



(24) 


where primes indicate differentiations with respect to x^ . 

From Fig. 4, we conclude that we must still account for the dis- 
tributed forces pyj and p 21 -. Regarding these forces as nonconservative, 
and assuming that they do not depend on the elastic deformations, we can 
account for their effect in the form of the nonconservative work 


nc 


[& ,♦ 

= (p . v. + p • w-)dx. 
Jq 71 1 v z-\ 


(25) 


so that the total work can be written as 

w = W c + w nc = -V + W nc ( 26 ) 

where the conservative work has been recognized as being equal to the 
negative of the potential energy. 

The system differential equations of motion, and the appropriate 
boundary conditions, can be obtained from the extended Hamilton's 
principle (see Ref. 1, Sec. 2.7) 

| 2 (ST + 5W)dt = 0 (27) 
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where all the virtual displacements must be set equal to zero at 
t = t«j stg- Introducing the Lagrangian L = T - V, Eq. (27) can be written 

as 

2 (6L + W nr )dt = 0 ( 2S ) 

J tl 


where the Lagrangian has the functional form 


L 




j = 1,2,3 ; i = 1,2 


. ,n 
(29) 


It will prove convenient to separate the Lagrangian into that 
associated with the rigid domain D 0 and those associated with the elastic 
domains . Hence, let the Lagrangian have the general functional form 
(see Ref. 19) 


L(t) « L 0 (t) + E 


1=1 


L.(x i ,t)dx i + L^-.t) 


(30) 


where 


Lq ( t ) = LoCQj ( t) » Qj (t )] , j = 1,2,3 


L 1 - (x^ ,t) - C 0 j (t) (t) ,v^ (x 1 . ,t) , v 1 -(x i> t),v](x i ,t),v^(x i ,t), 

w i (x i ,t),...,w^ , (x 1 ,t)] 

i = 1 ,2,. .. ,n 


(31) 


(32) 


L i ( # t) = L i [8j (t) ,63 (t) , v-j ,t) ,v i (£ i # t) ^ Ui ,t) f w i (^ ,t)] 


in which Lq is the Lagrangian corresponding to the system in undeformed 
state, Cj the Lagrangian density associated with any point of the elastic 
member i, and L^ the Lagrangian corresponding to the tip mass. Moreover, 
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jfc. represents the length of member i. From Eqs. (30), (31), and (32), we 
conclude that 


3 

6L = l 

j=l 


66. + |r- 

39j J 36j J 


n i 

r£ i r 

+ l - 

i=l 

Jo ^ 


aL- 


aL-; . aL- 

<5V-{ + ~Z7T~ <$V- 


6v- 


3L-j 


[a^7 ov i T a^- i avj 

at-i 


dx^ 


3L4 •) aLi ati * a 1 -! 1 , 1 in 

+ 37V 6V i + ‘3w7 5W i + 5w7 6w i + 3^" iW i 3WV 5w ij 

+ 6V ^ ’ t) + 3V^TTtT 6; i (a i ’ 4) + sw ^ 


3L^ 


+ 3w i (t. i ,t) ' iw i<*i -O 


(33) 


In addition. 


{W nc = }q ( Pyi 6v i + p zi 5w i )dx i 


(34) 


Inserting Eqs. (33) and (34) into (32), and integrating by parts with 
respect to t, we arrive at Lagrange's equations for the rotational m otion 


— - i!r f%~] = 0 • J = 1 - 2 > 3 

36j dt ^afijj 


(35) 


Moreover, integrating by parts with respect to t and x^- , we obtain 
Lagrange's equations for the transverse displacements , and the associated 
boundary conditions, in the form 


faL-j) 


sL.,- 


av^ at (avjj 


ax. 


till + £N + 


(3v J 


ax? UvTJ ry1 


0, 0 < x.j < , 

i = 1 ,2,. . . ,n 


(36a) 


and 
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3 

r 3L-j 

— -] - 

3vj 

0 

3 X-j 

NJ 


9V- 


3_ 

at 


3L. 


l 3V U 


aL ” 

sv i = °» ^ 5v i = 0 at x i = £ i 


i = 1,2,. ..,n 


(36b) 


3L i 

avT 


a 

3X: 


av H 


'll 


aL 


<$V: 


= 0, - 3^- «v! = ° at x,. = 0 


L uJ 


Equations similar in structure to Eqs. (36) can be written for w i by 
simply replacing v i by w^. 


3. Nontrivial Equilibrium 

Let us consider the case in which p yi = p zi = 0 and define an 
equilibrium configuration as a set of dependent variables e j , v 1 , w i 
constant in time and satisfying Lagrange's equations. Because these 
variables do not depend on time, they must satisfy the equations 


!L=0, j = 1,2,3 

36J 


and 


_3_ 
3Vj " 3X i 


faL^ 

T 


3 VJ 


2 (aL,! 


lav? 
i 


ax, v i 


= 0, 0 < x, < ij , i - 1,2, — ,n 


(37) 


(38a) 



[^1 

3L i" 

3V- ~ 3X 1 

3V- 

k. ' ^ 

3v.j 


3L i 


&v i = 0, ^TT 6v: =0 at X i = z. 


i = l,2,...,n 


(38b) 


__ 3 _ 
d\f\ ' 3X-j 


faL i 1 


dv!j 


^L i 

6v i = ~ dVf 5v i = 0 at x i * 0 


as well as a set of equations similar to (38) for w.. We shall denote the 
solutions of Eqs. (37) and (38), together with the set of equations for w. 
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by e. 0 , v iQ (x i ) > w.qCx^, where the first are constant and the latter 
functions of the spatial variables x i alone. . 

4. Perturbations About Equilibrium. The V ariational Equations of Motion 
The interest lies in the stability of the system in the neighbor- 
hood of the nontrivial solutions 0j O , v i0 (x.j), w io^ x i ^ * We sha11 S6ek 
stability criteria by means of Liapunov's direct method, and, to this end 
we let the solutions of Eqs. (35) and (36) and the companion equations to 

(36) have the form 


6j(t) = 0jQ + J “ 1*2,3 


V-j (x.j jt ) 0 ( X-j ) 


+ Vi1 (x i9 t), W 1 (x i ,t) 


= w i o ( x i ) + w il * 

i = 1 ,2,. . . ,n 


(39) 


where e^U), v^x^t), w.^x^t) are small perturbations. Inserting 
Eqs. (39) into Eq. (30), and expanding a Taylor's series about the non 
trivial equilibriums we obtain 


L - L ( 0 jo ,v iO ,v iO ,v iO 5 W iO >w iO* w io) + ^ 6 Jl 30 jo ^ 


n 

+ l 
i=l 


&i 


\\ . 9L • • 8L<? „ 3L-? 

9^ V il + 3^ Vi l + ^?0 Vl1 + ^ Wil 


W 1 + ... 


8L i 

+ w i'l 


dx^ + 


3L i 


3L i ■ 


3 V 


^ V il + 


X-j J 


1 3 3 { a 2 L 


A 


a 2 L 


+ 2 E z 


Vkl + 2 =lu Vkl + 8jl ° kl 


j=l k=l I 38 j0 3e k0 


J j0 kO 


jO kO ° i 
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r 3 2 L 


i .2 


3v iO 


" v ii + — ^ v 


3V 


10 


1 n 
+ i- £ * 

* 1=1 J 

l 

a 2 !., - 2 a 2 L i ,2 
+ ^ W " + ^ Wi1 ^ 


Ai -2 + Ai v ,2 + Ai v „ 2 + !!iL w 2 

“ V-, + V . 1 + —£ V il + i|<|2 W n 


IW Ml 

3V i0 


3V 


10 


3W i0 


S^L* 2 3 2 L. o ^ 

1 j- I w'.'T + 2 E 


’ll 


Ai 


j=l I 39 j0 3v i0 




Ai 


9 2 Lj 


3 2 U 


3 




+ ... + 


39 jO Sw 10 Jl 


f 3 2 U 


+ 2 


I 3v i0 3w i0 

3 2 L 


b v n"n + + ^T 0 ^ + "• 


a2 h 

'i0 3V iO 


a 2 U 


■ -j vUv'.’i + - w 4l W i 1 

av «' "1 1 ll »w!_swV_ 

3v i0 3v i0 


9w 10 9w i0 


dx^ + 


r 3 2 Li 2 3 2 L i - 2 
— O 1 V?1 + —f- Vil 

3v i0 10 


3 2 L, 2 a2l -i '2 3 

+ —T- w ii + w n 2 \ 

■2 1 1 3W| 0 1 J=1 


3W |0 


3 2 L- 


b 2 l . 


39 j0 3v iO 


e ji v n 


39 j0 3 Vi0 


e jl v il + 


3 2 Lj 


39 jO a ”iO 


3 2 Li 


e j1 w il 


+ 2 


3 2 L. 


3 2 L 1 . - -i 

V W + — * 7 V>-jW-i 

il 11 8V i0 3W i0 11 


I 3V i0 3W 10 


3 2 L^ 


'll 


+ J w • t v • i + — r — V.* -iWj -» 

3 w i0^i0 11 1 3v i0 9w i 0 1 


x. = iti 
1 1 


+ ... 


(40) 


where 3L/30.Q - 3L/30j 


e j * e j0 s 


, etc. But the term L(GjO’ v iO’* • • > w iV 


is constant. Moreover, by virtue of Eqs. (37) and (38) and the companion 
equations for w i , all the linear terms in the perturbed variables in 

expansion (40) reduce to 


aL 


j=l 3S j0 


Hi + E 
Jl i=i 


0 


aL. 

i 


l 9 Mo 


v n 


aL. 

i 

a^io 


w il 


dx. + 

i 


faL- 


^ Vil 


20 



3L . 

+ rr — W---J 

3w i0 "l x, « t. 


which are all linear in the generalized velocities ejpV^-j ,w^ . In view 
of this, if we retain terms through second order only, the Lagrangian 


becomes 


L = T 21 + T U + T 01 


where 


3 3 


= 2 E 1 • 0 jl 8 kl + 2 *• 

1 .1=1 k=1 38 j0 38 k0 1=1 


•S- p 2 h -2 + .2 

0 ^” Vil 9W 2 11 

0 _8V i0 SW i0 


3 2 L, • • 3 3 2 L 1 • • a 2 L 1 - • • 1 

+ 2 -= 7 — V..W.. + 2 E — 3 — 0..V-, + -7 — -7 — e il w il dx. 

3v i0 3w i0 1 11 j=l 30 j0 sv 10 J 11 3e j0 3w i0 3 1 J 


Vu o 3 2 L ■ 2 3 2 L,- • ■ 3 

+ ~ vf, + -Tp 2 -*., + 2 — + 2 E 

L 3v iO 3w i0 3v i0 3w i0 


3 3^L 1 


j=l |3ej 0 av i0 


e jl v il 


9^L.j 

+ a6 in^7n 6 J lW ' i ' 1 

JO i0 j x. = 


is quadratic in the generalized velocities 


2 2 3 2 L 


n (, ti r 3 r »<L. 


T n = E E — 2 — t- — 6 il®kl + z i z Ih — W~ 9 ^ v ,-! 

11 j=l k=1 3e jO 3S k0 J 1=1 Uo [j=l v 8 j0 5v i0 31 11 


S 2 ^ . 3 2 L- . 3 2 L i . 

38 j0 3v 10 J 30 jO 3 *iO J 11 3e j0 3w i0 


jl w il 
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3 2 L. 


9 2 L. 


— v n*n + - 


3Vj()3W-jO 


3 2 L- 

— e jl v il + 


3v i0 3w i0 

3 2 L- 


V il W il 


■n 

r 3 

dx^ + 

E 

>1 


a 2 L. 


e -t v 


36 j0 3v i0 


99 jO aw i 0 
3 2 U 


0j 1 W il + — 


3 2 U 


9 . v ’-■j 

+ — V.-.W.-. + — : v il w il 

3Vio3w i0 3Vjo^ w io 


30 jO 3w iO 


*1 = *1 


I 3e j0 3v i0 


0 jl w il 


jl il 


( 44 ) 


is linear in the generalized velocities, and 


3 3 3 2 L n 

T 01 ' V 1 = A ae jO 30ko e j l6kl + 1=1 





f ,f [- 

J 

l 

1 

o 


4 2 ^ !-• 2 
^ V il + 


3 2 Li ? A, 2 3 2| -i .2 3 ^ L i ,.2 „ ^ 

+ \ v »i + — 5- w., + — r? w.!, + —^2 W 11 + 2 . E , 

3^ 11 3W? 0 1 3w io 3w io j =1 


3 2 Li 


3 2 L, 


3 V w io 


— ejiw,! 


10 

/ A 

f 3 2 Lj 


l30jO 3v iO 


e jl V il 


3 2 L, 


+ 2 [i^i^ v ii w ii + 3v: 0 3vv o v n v n 


3 2 U 


’ll 


3w iO* w iO 


2 — w‘ w" 
il il 


r.2 


dx. + 

l 


3 £ L 


2 ^ a 2 l- 1 2 

1 v., + L w. n + 


3 f 3 2 L • 


v., + w., 2 E 

3v i0 11 0 11 ^ 


0 V. 


l 30 jO 3v iO 


jl il 


A, 

+ 0 .,W., 

30 jO 3w iO Jl 11 


X 1 = l i) 


(45) 


is free of generalized velocities. 

In view of the above, the perturbed Lagrangian can be written in 1 

the general functional term 


L = L ( 6 ji > 0 jl » v il ^il » v il »• • * ,w il * w il ) * J = ^ ,2 *^’ 


i = 1 , 2 ,. . . ,n 


(46) 
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Consequently, the variational equations can be written in the form of 
the Lagrange equations, Eqs. (35) and (36), but with the subscripts j and 
i replaced by jl and il, respectively. Unlike Eqs. (35) and (36), the 
variational equations possess trivial equilibrium. 

5. Discretization by a Ra.yleigh-Ritz Approach 

The variational equations discussed in the preceding section con- 
stitute a set of hybrid differential equations, in the sense that the 
equations for the rotational motion are ordinary differential equations 
and those for the elastic displacements are partial differential equations, 
where the latter are subject to given boundary conditions. It will prove 
convenient to transform the system into one consisting of ordinary 
differential equations alone. This can be done by using a discretization 
procedure based on the Rayleigh-Ritz approach. Indeed, let us introduce 
the notation 


Qj-j(t) - qj (t) , j - 1,2,3 


vn(xi »t ) 


p+3 

s (x-t )q 4 ( t ) , 
j=4 J J 


2p+3 

w.-tx-.t) * ■ X >l' i (x 1 )q i (t) 
11 1 j=p+4 J 3 


3p+3 

v,,(x,,t) = x <t , -i(xo)q • (t) , 
21 2 j=2p+4 J J 


4p+3 

w,i(x,,t) = X *i(x ? )q.(t) 
21 2 j=3p+4 3 2 J 


(47) 


(2n-l)p+3 

v,(x„,t) = x <f i (x n )q,(t), w , (x ,t) 
nl v n j=2(n-l )p+4 3 " 3 n1 n 


2np+3 

z 

j=(2n-l )p+4 


^(XpJqjft) 


where 4>-j(x.*) and (x^ ) are admissible functions, taken as the eigen- 

J i J I 

functions of the linearized system. With this notation, Eq. (43) can 
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be written in the matrix form 


T 21 = l {q(t)} T [m]{q(t)} 

where [m] is a constant symmetric matrix having the elements 

, j ,k = 1 ,2,3 


3 2 L 


ffljk = 3§ j0 3 \o 


m,,. = 


3 2 L, 


, 1 — ]— » k (x 1 )dx i + 

Jk J 0 30 jO 3v iO 


r 3 2 U 


r^7-*k< x i> 


L 30 jO av iO 


x i - 4 i 


(48) 


(49a) 


j = 1,2,3; k = 2(i-1 )p+4, 2(i-l )p+5, . . . , (21-1 )p+3, 

i ~ 1,2,... ,n 


m 


jk 


= -j — r* — ^ k ( x i) dx i + 

In BB in dW in k 1 1 


0 <* y jO dW iO 


r 3^L. 

-r i— ♦ k (x 1 ) 

39 j0 3W i0 


X. = Zj 
1 1 


(49b) 


j = 1,2,3; k = (2i-l)p+4, (2i-l )p+5,. . . ,2ip+3, i - l,2,...,n 


m jk = J 


f*i 3 2 L. 

~rj (x^ )*k( x i )** x 1 ) + 

0 3 Vl g J 


■ 

r i* j (x i )* k (x i ) 


9V 


10 


(49c) 


x i ~ M 


j,k = 2(i-l)p+4, 2(1-1 )p+5 , . . . ,(2i-l )p+3, i = 1,2,...,n (49d) 


.^i 3^L' 

m.. - ~ <j> j (x^ )'l'k(x-i )dx. + 

Ok J 0 3v i0 3w i0 J 11 


1 3 2 Li 

-t — r- — 4> i (x 1 )>i>k(x i ) 
3V i0 3W i0 J 1 K 1 


X i = *1 


j = 2(i-1)p+4, 2(1-1 )p+5, — ,(2i-l)p+3, 

i = 1 ,2 n (49e) 

k = (21-1 )p+4, (21 -1 )p+5,. . . ,2ip+3, 
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ra jk = 


r 1 ' A, 


ri +j( x i)*k (x i )dx i + 

iO 


q 9W^ - ^ 


3 2 L 


L 3W ! 


rsr*j (x j)+k (x i) 


x i = *1 


j,k - (21 -1 }p+4 , (2i-l )p+5,. . . ,2ip+3, i * l,2,...»n 


(49f) 


On the other hand, Eq. (44) leads to the matrix form 


T-j 1 - {q(t)} T [fKq(t)> 


(50) 


where [f] is a constant square matrix with the elements 


f ... = 


3 2 L 


jk ae j0 3 ^k0 


, j,k = 1,2,3 


r^i 


f jk 


r ,2 


0 39 j0 3v i0 


a a. 


L 30 jO 3v iO 


. •t' k ( x i ^ 


jX i 


= i. 


(51a) 


j = 1,2,3; k = 2(i-l)p+4, 2(i-l )p+5, . . . ,(21-1 )p+3, i « l,2,...,n 

(51b) 


f jk 


rSLi 3 ^L- 

I ] — ^(xijdx, + 


r 3 2 L.- 


— : — 'l>i c ( x i ) 

3e j0 3w i0 


x i “ *1 


j = 1,2,3; k = (2i-l)p+4, (21-1 )p+5, . . . ;2ip+3, 1 = l,2,...,n 

(51c) 


f jk = 


3 2 L. 


1 — — ♦,(x 1 )dx. + 

0 30 kO 3v 1O J 


r 3 2 L. 


2— ♦ i (x i ) 

30 kO 3V 1O J 


JXi 


j = 2(i-l)p+4, 2(i-l)p+5,...,(2i-I)p+3, i “ l,2,...,n; 

k = 1,2,3 


( 51 d) 
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f l 1 3 2 L} 

f ik = r? — r )dx- 

Jfl 3 1 1 




—7 ^ . (*,• ) 

ae k0 3w i0 J 1 


J *i = h 


0 d0 kO dw iO 

j = (2i-l)p+4, (2i-l )p+5,. . . ,2ip+3, i = l,2,...,n; k - 1,2,3 ( 51 e ) 


f jk = 


r£ i 3 L, 


r a^L.- 


. J — ♦i(x i )*lc(xi)dx 1 - + 77 — 4>i (X - )<» k (Xi) 

o 3v io aw i0 J _ 3v io 3w iO 

j = 2(i-l)p+4, 2(i-l )p+5,.. . ,(2i-l)r+3, 
k = (2i-l)p+4, (2i-l)p+5,, . . ,2ip+3, 


aX i 


= u 


i = 1,2,... ,n 


(51 f) 


: 3^U I 3-L_- 

f ik = I 4. k (x i )i(i . (x. )dx^ + — *|.(x. )<p i (x i ) 

JK 3v io aw io 1 J 1 1 jsv i0 3w i0 k 1 J 1 


r 3 2 L, 


x i 


- JL 


j = (2i-1 )p+4, (2i-l)p+5,...,2ip+3, 


k = 2(i-l )p+4, 2(i-l)p+5,...,(2i-l)p+3, 


i = 1 ,2,. . . ,n 


Finally, from Eq. (45), we can write 


T oi - v i = - ^q(t)} T [k]fp(t)} 


where [k] is a constant symmetric matrix with the elements 


3 2 L 


k jk “ 30^30 


jO kO 


, j »k = 1,2,3 


k jk 3 


■-! 


n-i S z L j 


0 30 jo 3V io 


-*k (x i )dx i ■ 


3 2 L - 


L 30 jo 3v io 


— *k< x i) 


x i - *i 


j = 1,2,3; k = 2(i-l )p+4, 2(i-1 )p+5,. . . , (2i-l )p+3. 


<51g) 


(52) 


(53a) 


i = 1 ,2,. . . ,n 


(53b) 
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*1 3 2 L 1 - 

0 3e j0 3w i0 ** 


r 32l i 

(x i )dx i ' aijoaiTjo' ’M x O ^ = £ _ 


j = 1,2,3; k = (2i-l)p+4, (21-1 )p+5, — ,2ip+3, i = l,2,...,n 

(53c) 


A * 

2 2 3^L 

1 MVMV + — ^ ♦j« x i ) *kf x i> + 7^i ^( x i^k^ x i^ 
Jo I 3 v2 q J 1 K 1 3v] 2 J 1 avVg J 


3 2 L 1 
L j. ii 


r k-k-ktx.,.) + ♦”(x i )*^(x 1 )j cjx. 


3v i0 3v i0 


2 ^ x i ( x i ) 

3v i0 J x i s *1 


j ,k = 2(1-1 )p+4, 2 ( i -1 )p+5 (21 -1 ) p+3 , i = 1,2,. ...n (53d) 


M 3Z h 

0 3v i0 3w iO 


■ g2[_ . 

-*j(x 1 .)* k (x i )dx i - av.V^::' *J< x 1>*k< x 1> 
iO i0 i0 J x = 


j = 2(1-1 )p+4, 2(1-1 )p+5,. . . ,(2i-l )p+3. 


k = (21-1 )p+4, (21-1 )p+5, — ,2ip+3 


i = 1 ,2,. . . ,n (53e) 


|3 2 L, 3 2 L i a 2 L, 

] — 7 * i (x-)* k (x 1 ) + —4 *;(x i )* k (x i ) + — J +jf x 1 )*k(x--i) 

0 3w,£ J 1 K 1 3w]£ 3 1 3wljn J 


3^L - 

+ aw' aw" fakkk> + K 

iO iO L J 

j 
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8 Z L- 


'J> w iO 


'Pj (x i )ip k (x 1 ) 


jX i = £ i 


j,k = (2i-l )p+4, (2i-l)p+5,...,2ip+3, i = 1,2 n 


(53f) 


Introducing Eqs. (48), (50), and (52) into Eq. (42), we can write 
the Lagrangian in the matrix form 

L - |{q} T [m]{q} + {q} T [f]{ql - ^q> T [k]tq) 

Using the approach of Ref. 21 (see Sec. 3-4) , we can write Lagrange's 
equations in the matrix form 


(54) 


d 

dt \aq 


© - © - <o » 


(55) 


Hence inserting Eq. (54) into (55), we obtain the equations of motion 

[m]{q> +([f] T - [f]){q) + [k]{q> = {0} (56) 

so that, introducing the notation 

[g] = [f] T - [f] (57) 

T 

where [g] is a skew-symmetric matrix, [g] = -[gj* we obtain 

[m]{q} + [g]{q) + [k]{q} - {0} (58) 

where [m] is identified as the inertia matrix, [g] is a "gyroscopic" 
matrix and [k] is a stiffness matrix which includes terms due to elastic, 
gravitational, and centrifugal effects. 
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6. Liapunov Stability Analysis 

We shall seek criteria for the stability of motion in the neighbor- 
hood of the nontrivial equilibrium by means of the Liapunov direct method. 
This is equivalent to the problem of stability of the perturbed motion 
about the trivial solution. In terms of the discretized system, the 
perturbed motion is described by the vector {q(t)}, so that the interest 
lies in a stability analysis about the trivial equilibrium {q} = (0). 

It was shown in Ref. 75 that the Hamiltonian is a suitable 
Liapunov function for the type of problem at hand. Assuming that the 
system possesses a certain amount of internal damping, however small, the 
equilibrium is asymptotically stable if the Hamiltonian is positive 
definite. In terms of the perturbed variables, the Hamiltonian has the 
form 

H * T 21 - T q1 + V 1 = l{q(t)} T [m]{q(t)}+ ^-{q(t) } T [k]{q(t ) > (59) 

But the function T 21 is positive definite in the generalized velocities 
q j ( t ) by definition. Hence, if the function 

•i 

k = f{q(t)} T [k]{q(t)} (60) 

is positive definite in the generalized coordinates q . (t } a then the 

J 

Hamiltonian is a positive function in the generalized coordinates and 
velocities and the equilibrium is asymptotically stable. The function 
k is positive definite if the matrix [k] is positive definite. Whether 
[k] is positive definite or not can be ascertained by means of Sylvester’s 
criterion (Ref. 1, Sec. 6.7). The matrix [k] will be referred to as a 
Hessian matrix. 
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7. Natural Frequencies of the Complete Structure 

The Liapunov direct method provides qualitative information con- 
cerning the stability or lack of stability of an equilibrium configuration. 
Similar information can be extracted from the system of equations (58) 
via the eigenvalues. In addition, the eigenvalue problem yields results 
of a more quantitative nature in the form of the system natural frequen- 
cies and the normal modes for the complete structure, where the latter are 
defined later. It turns out that Eqs. (58) lead to an eigenvalue problem 
of a special nature. The nature of the eigenvalue problem can be con- 
veniently discussed by converting the set of equations from second order 
to first order. Indeed, if the configuration vector (q(t)} is of 
dimension N» then we can introduce the 2N-dimensional state vector 
{x( t ) > in the form 


{x(t)> = 


r {q(t)}| 


(61) 


No confusion should arise from denoting the state vector by {x(t)>, because 
the symbol x^ used to denote the position of a point in the elastic 
members represents a spatial coordinate independent of time and not a 
time-dependent generalized coordinate. Accordingly, if we introduce the 
2N x 2N matrices 


[M] 


TM [oj 

10] [k]. 


' [g] IkT 

_-[k] to].. 


( 62 ) 


then the set of N equations (58) can be transformed into a set of 2N 
first-order equations having the matrix form 
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[M]{x(t)} + [G]{x(t)} = {0> 


(63) 


where [M] is symmetric and [G] is skew-symmetric, 

CM] = [M] T , [G] = -[G] T (64) 

because [m] and [k] are symmetric and [g] is skew- symmetric. 

The matrix equation (63) is of the special form treated in Ref. 20, 
so that the eigenvalue problem can be solved by the method developed 
there. Hence, letting 

{ x ( t ) } = e xt {x> (65) 

where X and {x} are constant, we obtain the eigenvalue problem 

x[M]{x) + [G]{x> = {0} (66) 

It is shown in Ref. 20 that the solution of the eigenvalue problem (66) 
consists of 2N eigenvalues x r and eigenvectors {x} f (r * 1,2,...,2N), 
where the eigenvalues consist of pairs of pure imaginary complex conjugates, 
x r = ± iw r , and the eigenvectors also consist of pairs of associated com- 
plex conjugates (x) r and {x*} r (r = 1,2,...,N). Moreover, the eigenvectors 
are orthogonal in a certain sense. Reference 20 provides an algorithm where- 
by the eigenvalue problem can be solved in terms of real quantities. The 
method will be used later in this work to solve the eigenvalue problem for 
a specific spacecraft. 

8. Lagrange’s Equations in Explicit Form 

Lagrange’s equations, Eqs. (35) and (36), are written in a 
general form. Before obtaining the nontrivial equilibrium and the 
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corresponding variational equations, we must express them in a form in 
which the various coordinates appear explicitly. By virtue of the assump- 
tion that the satellite mass center moves in a circular orbit with 
orbital velocity ft, we can replace K/R c by ft in Eq. (5). Moreover, the 
first terms in Eqs. (4) and (5) can be ignored because they are constant. 
In view of this, if we recall that the Lagrangian can be written as 
L ~ T - Vq - Vjr^ - V^ B , then we can substitute Eqs. (4), (5), (15), and 
(23) into L, and obtain 

L(t) = + {w}^"{K} + T E +1 ft^tr[J^)] 


- f n 2 U a } T [j(°>]U a } - 1 S f [P x1 (v! 2 + wJ Z ) 


+ EI zi vf(l - fv,: 2 ) + EI yi wV 2 (l - f wj 2 )]dx. 
where P x1 - is the axial force at any point of the slender rod, and 


(67) 


[J (0) ] « " [j| 0> ] = " 

i=0 i=0 


(68a) 


(K> = ^ Pl [h| 0) + r| 0) + + mi [hj 0) + r\ 0) 




x i = *i 


(68b) 


1 " 
TV = 1 i 
E 2 • i 
c 1=1 


f S/ * T 9 V 

1 p4{u.} {u }dx. + m.{u.} T {u.}| 

JO 1 1 ’ 11 1 •x. = i 


V 


(68c) 
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in which [J^)] is the inertia matrix of the body in deformed state in 
terms of the reference system xyz, {K} is an angular momentum matrix due 
to the elastic velocities, and T E is the kinetic energy due to the elastic 
velocities. The elements of [J^] are given by Eqs. (9) and (10). Intro- 
ducing the notation 


°ill = < 3 i[^yi + v 1 0 2+ (h zi *+w 1 ) 2 ] , J in ( ) “ m it( h yi +v i ) 2+ ( h zi +w i ^-1 


x i =£ i 


3,22 = P 1 -C{h xi +x i ) 2 +{h zi +w. ) 2 1 , J i22 (^) = ■ 1 C(h x1 +x,) Z +(h z1 « 1 ) 2 ]| 


0,33 - PiC(h xi + *i) 2+ (h y,+v,) 2 ]. J i33 (M) = -,i:<h x ,+x,.) 2 +(h y ,*v,) 1 


X .=£ . 
1 1 


J il2 " ^i 21 = " p i ^ h xi +x i ^ h yi +v i ) » = J i21^i^ ” “ m i^ h xi 

+ X, )( h y, + V, ) jx-sj,. 


J il3 ” J i 31 " “ p i ^xi +x i Hh zi -+ W 1 - ) , 3^i ^ " J i31^^ ~ ~ m i ^xi 

+x i)(h z i+ w i) 


O-j' 23 ~ ^i 32 = ~ p i< h yi +v i )<h zi +w i^ J i23^i^ * J i32^i^ = “ m i ^ h yi 


+v^ )(h zl *+w.) 


X-! £.4 

1 i = 1 ,2,. . . ,n 


(69) 


we can write 


[ 0 ,] 



[J,]dx. + [J,(t,)] 


(70) 
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In a similar way, from the second of Eqs. (68), we have 


{K} = l 1 {K. }dx. + IKAU)) 
i=l [JO 1 1 


( 71 ) 


In view of the above, the Lagrangian densities can be written as follows 


L 1 -<x i ,t) = 1 {u 1 } T [j!°h{u)} + {<u} T {K i 1 + 1 p 1 -{ii i } T {u i } + j ft 2 tr[5!°^] 


1 Jl, 


(0) 


- f ^{^ a > T C3S 0) ]U a } - l P xi (v-. 2 + w! 2 ) - 1 EI zi vf (1 - f v! 2 ) 


- i 'EI yi wV 2 (l - |w! 2 ), i = 1 , 2 ,. . . ,n (72) 

whereas the parts of the Lagrangian associated with the discrete masses 
are 

LiUj.t) = 1 {^[jj 0 ’^-)]^} + {(o} T {K.(^.)} + 1 m- {u-> T {u i > 

£ c x i =J ^i 

+ 1 S2^tr[j|° ) (t i )] - | JJ 2 {£ a } T [j|° ) (£ i )]{£ a } (73) 


From the context it should be obvious when brackets and braces denote 
matrices in Eqs. (67) through (73) and when they do not. 

Substituting Eq. (67) into Lagrange's equations for the rota- 
tional motion, Eqs. (65), we obtain 




[J (0) ] U>1 + 




- ft (jfr W T J ([J (0) ] M + {K}) = 0 j = 1,2,3 ( 74) 
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Moreover, Lagrange's equations for the transverse displacements Vj are 


i (aV7 [ b { ° ,] ) { “ } + {U)}T (air CK,»] + tr (^7 [ 3 | 0) ]J 


.1/ 

2 


U a ) T f_L. [j| 0 ) ]){Ji a } - -i [{U)} T _|_ {K} + Pi V,- 

a 3v,- 1 d at a vj ‘ 

V. | V I 


~ — ^xi v i + % ^zi v i v i 2 ^~^zi v i ^ 

2v . c. a v • 


~ v^ )J + Py^ = 0 , 0 < < £-f , i l,2,...,n 


which are subject to the boundary conditions 


■ (p xi - | EI zi b' 2 ) v i + afr CEi zi V? (1 - f vj 2 )] 

+ 1« 2{ W T 3^[4 0) ] { “ } + trg7T[0 {0) ] -3U a } T ^[ j(0) ]{V> 


- if «-> T air {K} + p i h) = 


EI zi v i U - I *J 2 > - 0 


at x.=A. f i 


v. - 0, vl = 0 at x 1 = 0, i - 

Similarly, Lagrange's equations for the displacements w. are 

? {W}T (a^ c3 i 0)] ] { " } + { “ }T (afr^1>] +\* (a^ c3 !° )] 
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|sj 2 U a } T 


_3 

3W 


r [i{ ol i 



r 

{w} —7— {K} + p-iWj 

3 w i 1 




w] 


+ I El . w'.wl' 2 )+ — 2 — g-C-EI - w. (1 
2 yi 1 1 sx,-^ J 1 1 


c 1 2 

- | w. )] + P Z1 - - 0 , 0 < x- < tf , i = 1 ,2,... ,n 


(76a) 


which are subject to the boundary conditions 


-( p xi - f EI zi «i' 2 ) w i + £7 t EI yi - I w i 2 )3 


+ l a 2 
2 


M J £ T [0S O) ] ♦ tr » [J< 0) 3 -3 U a J T ^ [J <0, 3t^> 


_ i- 
3t 


{u } T 4- {K} + Pi w, 

3W-j 1 1 


= 0 


Ely-j wV(l - | w! 2 ) = 0 


p at x . = & .j , i - 1,2,... ,n 


(76b) 


w- s 0, wl = 0 at x. a 0, i s l,2,...,n (76c) 

9. Equilibrium Equations in Explicit Form 

For a gravity- gradient stabilized satellite, the angles Oj(j = 
1,2,3) are measured relative to an orbiting system of axes. The orbit 
being circular, with the orbital angular velocity being equal to fl, the 
orbital axes rotate relative to an inertial space with angular velocity a 
about an axis normal to the orbital plane. This axis is denoted by c (see 
complete definition later). Hence, the angular velocity matrix {u>} can 
be written as 
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( 77 ) 


(o)} = + {w) r 

where U c ) = { Ji c ( ej ) > is the matrix of direction cosines between axis c 
and the reference system xyz, and (w} r = {o>( a j ,6j) > r is a matrix whose 
elements are the angular velocity components of system xyz relative to 
the orbital axes. They are linear combinations of the velocities Gj 
(j = 1,2,3). 

The equilibrium equations can be obtained by deleting from Eqs. 
(74) - (76) all the terms involving derivatives with respect to time. 

This implies that we can replace {u>} by in these equations. Hence, 

the nontrivial equilibrium must satisfy the general equations for the 
rotational motion 

« C } T [J (0) ] -f- {* c > - 3U a } T [J (0) ] t§- u a } = 0, j = 1,2,3 (78) 

90 j j 

as well as the boundary-value problems defined by the differential equa- 
tions 


\ « 2 [{*C> T !V7 Co5 0) 3t£ c > + tr ^7 CjS 0) ] - 3{2 a } T [jj 0) ]{£ a > 


+ - 9 - r(p • - - ei . v'. ,2 ivn — ^ fei -v-d - - V ! fc )i * o 

3 X - u XI 2 Zi v i * J ax 2 L zi 2 1 


5 ..,2 


0 < x ^ < 5,*, i ~ 1,2,... ,n 


(79a) 


and the boundary conditions 


-( p xi - I EI zi v i 2 > v i + 1x7 [EI zi v i' (1 - I v i 2 n + \ n 2 


ty T 7V7 t J i 0) ^* c > 
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1 


' at x.=£|, i = 1 ,2,. .. ,n 

(79b) 

Vj s 0, v| = 0 at = 0j i = (79c) 

Moreover, it must satisfy a set of equations similar in structure to Eqs. 
(79), but with v.j replaced by w. . 

10. The Variational Equations for the Discretized System 

The variational equations for the discretized system were obtained 
earlier in the form (58), where the matrices [m], [g], and [k] are defined 
by Eqs. (49), (51), (53), and (57). Although the equations just mentioned 
have the advantage of revealing the symmetry of [m] and [k] and the skew- 
symmetry of [g], the formulas for deriving the elements of the matrices 
are not the most suitable from a computational point of view. Indeed, we 
wish to present a procedure whereby the actual derivation of the variational 
equations is performed by a digital computer. 

Consistent with earlier notation, we shall denote quantities 
associated with equilibrium by the subscript 0 and perturbed quantities 
by the subscript ]. With this in mind, we can write the Lagrangian in 
the form 

L = L 0 + (80) 

L 0 = J M J 0 t j(0) lo M 0 + M r o<Vo + T E0+I n 2 tr [ j(0) ] 0 

- f {i a }J[J (0) ] 0 U a } 0 - V E0 .(81) 


+ tr^ 7 [ J (0) ]-3U a } T 3^[jJ° ) ](V 


(0) 


= 0 


5 „,2> _ 


- EI zi v lO -fvp = 


0 
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and 


= {w}|[J^]q{w}q + J {u}q[J^^]-j {w}q 

+ + m}{K) 0 + {u}{{K}-| + {w}J{K> 1 + T e1 

+ I n 2 tr[J (0) ], - 3S5 2 U a }J[J (0) ] 0 U a } 0 - f o 2 {» a }f[J (0) ] 0 U a } 1 


- - 3n 2 tJt a >|co (0) a 1 {A a } 0 - V E1 


(82) 


in which 


{to} " {o)}q + {u>> 1 = fK£ c > 0 + Z 


r f ) [ 
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"TO 


U) 


«1 + 


"iO 


{w} 


1 


3 3 


+ ~ 2 l 

2 1=1 j=0 I_| > 9e i0 3e j0 


{to} 


Wj + 2 


[ 3e jO 30 iO 


-{u} 


q j q i 


(83) 
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{£ a } 0 + s 
a U 1=1 


36 i0 




3 3 | „2 


+ \ l z Lq - K g — u a >l 
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Cm] 
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v il 
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(85) 


n frij „ 

m = z {K.-JdXi + {K,- (*.,-)> 
i=l [Jo 1 


{K} 0 = 0 


{K}, = 


- z P 

1=1 JO 


IV,; 



r* 1/ 1 


tf.o + 


11 T 1 


3K 


3W 


10 


w il + 


faK i 


l 3V 10i 


v il + 


3^. 


3w io 


w ii 


+ 


2 * 
JIj 


3v 10 3w 10 


> v il w il + < 


3?K. 


3v io 3w io 


v il w il 


dx^ + 


iO 


,3K.(£.)i 

Hn — } v ii^i) 


+ 1 


9K. (£. ) 

3 XI 


3W 


iO 


K-l ( *i> + ilTT- 


10 


K<*i) + I 




3w i0 


W^(llj) 


f3 2 K i (£ j ) 


3v 10 3w i0 


} V 11 (i 1 )w! 1 (t i ) + 


f 3 2 K i (£ i ) 


3v i0 3w 10| 


v i 1 ( £ i ) w 1l (J^l ) (86) 
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Note that Eqs. (83) - (86) represent Taylor’s series expansions of {<u}, 
U a >, [J^], and {K} about equilibrium. Moreover, we have 



1 

2 




( 88 ) 


as well as 
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To obtain the variational equations in terms of the discrete co- 
ordinates q j ( t ) {j - 1 ,2,. . . ,2np+3) , we must insert the modal expansions 
(47) into L-j and perform the indicated integrations over the spatial 
variables (i = l,2,...,n). Because the resulting expressions are very 
lengthy, we shall not write them explicity, but proceed with the derivation 
of Lagrange's equations instead. To this end, it will prove convenient to 
denote constant terms by the subscript c and terms that are linear in the 
generalized coordinates q^(t) and generalized velocities q.(t) by the sub- 

j j 

script a. This enables us to write 


!Lu 

aqj 


3 , ,T 

IT- 


h t 


[j (0) ] 0 (“} 0 
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(91a) 
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3Ln f 
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j = 4,5,...2np+3 ( 91 d) 

which enables us to write Lagrange’s equations for the perturbed motion 
in the compact form 
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As mentioned already, the advantage of Lagrange's equations (92) over 
those derived in Sec. 4 is that Eqs. (92) permit automatic derivation by 
means of a digital computer. 

Before specializing the equations to a particular satellite, let 
us derive an expression for the axial force in terms of matrix nota- 
tion, The axial force P X1 - is due to centrifugal and differential gravity 
effects. Introducing the modified potential energy density associated 
with member i 

V-l mod = - I + ‘ 3U a } T [jj°^]U a }) 

- 1 a 2 ({£ c } T [jj°hx i )]U c ) + tr[jj°^(x.j )] 

- 3{£ a } T [jj° ) (x i )]{^ a })6(x i - A, ) (93) 


where the terms inside parentheses and multiplying <$(x.j-£.j) are due to the 
tip masses, the axial force density can be written in the form 


xi 


(x.j) = - 


3V 


i mod 
9x,* 


= l a ZUlri T tii 0 ) VU r } 


tr[j| 0) ]‘ - 3{a a > T [5| 0) ]'{a a }) 


(94) 


in which we introduced the notation 

(o)-, . 3 r i(0) 


=^[3'°h + ^[3 i (0) (x.)]6(x r , i ) 


(95) 


j 2 , o 

Observing from Eq. (67) that P X1 - is multiplied by (v^ + ), we ignore 

any transverse terms in [J 0) ]‘, SO that using the first of Eqs. (68) and 
Eqs. (69) we obtain the approximation 
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Inserting Eqs. (95) and (96) into Eq. (94), we can write the axial 
P xl * at any point x- in the form of the integral 
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It follows that the desired expression has the form 
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1 ] . The RAE/B Satellite. General Formulation , 
a. Equations of motion 

Next let us specialize the equations to the case of a satellite 
consisting of a rigid core with six flexible booms, as shown in Fig. 6 . 
First, we wish to determine the matrices [£j] of the direction cosines 


between axes x.y. 

ri 

z. and 

xyz. From Fig. 6, it 

is easy 

to 
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Ooo) 


where sa = sin a, ca = cos a, S3 = sin 0 , and cb = cos 3 . Moreover, to 
write the angular velocity matrix {w} in explicit form, we must specify 
the rotations 8 j (j * 1,2,3). Assuming that system xyz is obtained from 
system abc by means of the rotations about y, -e-j about x, and 03 
about z, and recalling that axes abc rotate about c with the constant 
angular velocity ft , matrix {w} can be shown to have the expression 
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where se*j = sin e-j , ce-j = cos e-j s etc. Because the direction of the 
radius vector R c coincides with that of axis a at all times, the direction 
matrix {£ g } can be written as 


{£ a l 
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C0£ C03 - S0-| S02 S03 

- -(c0 3 S03 + S0] S02 C03) 
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(102) 


It will prove convenient to rewrite matrices {w} and {£ a > as follows 
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s&i 

C 0 -J 

. 


L C 0 2 ] 


(104) 


Introducing Eq. (103) into (4), and recalling Eq. (68), the kinetic energy 
becomes 
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T = ^ {u)} T [0 (O) ] {<4 + {oo} T {K} + T e 

= \ {0} T te*]{ [0]| fJ (0) ] [e] 3 [e*] 1 {83 + {e) T [e*]{ [e]] [J (0) ] {9> 3 

+ \ {S}J [o {0) ] {0} 3 + (9} t Ce*3 3 [033 00 + { 6 }J {K} + T E 

+ 0 {0}J [ 0 ]{ [e]J [J<°>] [ 0 ] 3 [ 0 *] 1 {6} + n { 0 }J [ 0 ]{ [ 0 ]J [J (0) ] {0} 3 

+ Q {e}J [e]{ [e ] 3 {K} + ^- r 2 {e}£ [e]j [e ] 3 [J^] [e ] 3 [e ] 1 {e } 2 

= \ {0) T [e*]j [J*] [0 *]t m + {0} T [s*]] [e]| [J] {e} 3 

+ {0) T [e*]{ {K*} + 1 J 33 e 2 + K 3 0 3 + T e 

+ a {e}£ [e]J [J*] [e *] 1 { 0 } + a {e ) 2 [e]J [e ] 3 [J] {e ) 3 

+ a {fl}| [ 0 ][ {K*} + la 2 ' { 0 }J [e]J [0*]. [e ] 1 {e } 2 (105) 

where [J*] = [e ] 3 [J^] [e 3 ] and {K*} = [e]I (K). Moreover, inserting 
Eq. (103) into (5), and recognizing that U a }^ = - ({0> 2 )^ [e]T (9)1, 
we obtain the gravitational potential energy in the form 

V G = - 1. R 2 tr [J (0) ] + | R 2 (-{0}^) T [0]j [9]J [j(°>] [e] 3 [e], {-(8} 2 ) 

= - 1 R 2 tr [J (0) ] + | R 2 ({0}^) T [e]j [J*] [0], {0>£ (106) 

where primes indicate differentiation with respect to eg. Expanding the 
matrix involved in T and Vq, and recalling that L = T - Vq - - Vgg, the 

Lagrangian L can be written in the form 

1 = 1 [J*-j e| + (J^ c^e-j + J ^3 s 2 e-j + 2 0^3 s 0 i ce i) ^ 
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,* :2 


- 2 (Ji2 c0 l + J 1 3 S0 1 J ®1®2 " 2J 13 0 1 0 3 + J 33 e 3 
+ 2 ( J 23 C0 1 + ^a 59 ]) 9 2 9 3 ] 


- K*e-| + (K^cBj + K^se^ e 2 + K 3 e 3 + - ^ (uT + vt + wpdm i 

+ q . {( J -| iS 02 + J -| 2 S0 - JC62 - J - J 3 c 9 ice 2 ) e -| + L - J ^ c 0 -j S 02 

- d * 2 S 0 1 ce -| ce 2 + J * 3 56 ^ 06^02 - J * 3 se^seg + 0^3 c0 2 

- s 2 0-j)]02 + {O33 C0-JC02 - 0*3 S02 - O23 se l ce 2 ) ®3 " K * se 2 

- <2 se - jC02 + K3 ceice2 } + ft 2 ( 0 ^ + O22 + J33 ) + Jr J]-j ( s ^ e 2 

- 3c 2 02) + j O 22 s 2 ®! ( c ^ e 2 “ 3 s 2 02) + 4 0* 2 S0-|S02C02 

+ 1 J * 3 c 2 e 1 (c 2 02 " 3 s 2 e 2 ) - 4 J *3 c 0 -|S 62 c9 2 + 0^3 se^ce] (3s 2 e 2 


.2 .2 .2 


t - ■» ^ r£-j 

e 2) ] - l e J P xi 


- C 


3 V. 


. I 2 


U 3x 'i 


* fc . 

“i 


^ 2 i 


dx ^ 


1 6 l Zi 

? 1=1 Jo 


El • 
zi 


fa 2 v i l 2 r 


l 8x i J 


1 - 


fdvj ) 


l 8x i 


2 -t 


+ El. 


fz 2 * j ) 2 


yi 


i 3x i 


3 W . 


t 3x ’j 


>dx.. 


(107) 


To obtain Lagrange’s equations for the rotational motion, we intro* 
duce Eq. (107) into Eqs. (35), and obtain 

[(J 33 - J 22 ) s 0 l ce l + J 23 e l " s 8 l ^ e 2 + ( J 12 S0 1 * J 13 C0-})0-|6 2 
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+ (-J23 se 1 + Jj 3 08^6263 - (K* SB] - K3 ce-j)e 2 
+ < J 12 ce*jC02 + 0*3 s e-j C02 ) 6] + [0*2 s0 i se 2 - ^*3 ce-|S02 

/ Vkr ^ O a 1^ ( 

+ (O33 - J22 ^ f C ^ 0 1 ~ 5^7)003 - 4 J23 ce|se*jc 92 ]S 2 

- (J33 S 0 ]ce 2 + O23 C 6 ]C 02)03 - ce*|C 02 - K3 S 9 ] C 02 > 

+ fl 2 [4 J* 2 CO] S02C02 + ( J ^2 " J33) S0]C0| (c^02 - 3s 2 02) 

+ ^ ^13 s9 1 s 0 2 c6 2 + ^23 ^ c ^ 9 1 - s^ 6 -j)( 3 s ^02 - ) ] 

- CJ *1 - < J *2 C0 1 + J 13 s0 l ) 0 2 “ ^*3 0 3 - K * 

+ n (0*7 se 2 + J* 2 S0ice2 - 0*3 007002)] = 0 ( 108 a) 

{[J*l 002 - J*2 se l se 2 + J *3 c 07 S 02 ]e 7 + [-0* 2 ce ] c0 2 + J 22 se l s0 2 c0 l 

* J33 C87S87S02 - J73 $87062 - J23 S0 2^ C ^ 0 1 ” s ^ e l)]®2 “ (^33 c0 7 s0 2 
+ J73 C02 - J23 s 8 *j s 0 2 ) ^ 3 “ K* c0 2 + ^2 S0 1 S0 2 ~ ^3 

+ ft 2 [4 0*7 s©2 C 02 + 4 0*2 s67(c 2 62 - s^e?) - 4 s 2 07 S 02 C 02 

- 4 O73 c07(c 2 02 - s ^ 6 2 ) + 8 j| 3 s67ce7S02ce2 - 4 O33 0^7 562062] 

" 3 F ^ J 22 c 20 7 + O 33 5% + 2 J 23 s e-j ce-j 3©2 - (J *2 C0 1 + J *3 s 0 l) 0 l 

+ (O23 ce l + J 33 se l ) 0 3 + Kj C07 + K3 $67 + fl [- 0 * 2 067563 

- O22 S 07 C 07 C 02 ** O73 S07S02 + (0^67 - 5^07)063 

+ 0^3 507007002]} =0 ' ( 108 b) 
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- J *2 ®1 + ^ J *2 C ^ 6 1 + J *3 se ] C0 1 )®2 + ^ J 22 “ J *l) ce l + ^23 s0 1^1®2 

+ O23 e-j e 3 + J* 3 00^63 + K£ 9-j + K* C0-|0 2 
+ to {[-2 J*2 se 2 + (J*-| - J^) se l ce 2 + J 23 c0 1 c0 2^1 
+ C(0?2 - J*-j) c0‘jS02 - 2J*2 se-jce-jcej + J23 S0 1 S0 2 
+ J*3 (c 2 e-| ** s 2 ©i )ce 2 D ©2 + C J 23 se 2 “ J*3-seice 2 )&3 + S02 

- K-j s0.iC0 2 } “*■ to?" , t~^|2 ( s ^ 9 2 ~ 3 c^ 02) 4(J-j-j - ^22^^1^^2^ 0 2 

+ 0*2 S ^ 0 1 (c 20 2 “ 3 s^ 02) + 4^23 c0ise 2 ce2 

- J*3 s9-jC6 1 (c 2 02 - 3s 2 62)] - ^ [-0*3 6] + (J* 3 C0] + O33 se-j)^ 

+ J33 63 + K 3 + to {-0*3 S02 - J23 s0 1 c0 2 + J 33 c0 1 c0 2)^ " 0 (108c) 


Considering Eq. (36) in conjunction with Eqs. (69), (99), (100), 
and (102), and letting i = 1 , we obtain the differential equation for v-j 

p 

p] <( h xl + X 1 ) [(w]Ca + u) 2 Sa) (a)] Sa - u^Ca) - 3u (A^ca + ^ a 2 Sa ) (^al Sct 

- JJ a 2 Ca)] + (hy-j + v-j ) [(oj-jCa + a> 2 Sa)^ + go 2 + 2S? - 3ft 2 (^al Ca 

? 2 \ 2 
+ i^sar + i a 3y] + (h z -j + W-{ ) [013(01-1 Sot - 0> 2 Ca) - 3ft ^a3 ( A a l s « 

- ^ a 2 Ca )] + w ] ( Cata>i + S aai2 ) } - [- (w-j + h 2 j)(u^Ca + W2 S °0 

+ (h xl + X-j) w 3 + 


3 

3X-| 


< 7 7 2 

- to {(£ c ”jSot - £ c 2 Cot ) + *C3 + 2 “ 3[(^ a -jSa - S' a 2 Ca ) 
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+ x * 1 " + l l ' ! h xl + + m ] (h x ) + ^l)J v| 


+ ? EI zl v 


W 2 } + z! [ - EI zi v ‘i ( ] - f v i 2)] + Pyl ■ 0 

Oat 


O < X-I < Z] 


(109a) 


which is subject to the boundary conditions 


v-| (0) = 0 , v}(0) = 0 (109b) 

o 

C(h x , + &]) [(w]Ca + a> 2 Sa) (a)jSa - wgCa) - 3 « (A a -|Ca 

I 

+ fc^Sa) (^ a -|Sa - 5 , a 2 Ca )] + (hy-J + V-j) [(wj Cct + u^Soi)^ 

+ W3 + 2 _ 3 ^2 ^£ a ^ Ca + ^ a 2Sa)^ + ^93^ ] + (h z i + W-j ) [0)3(01] sa 

- a> 2 ca) - 3 O 2 * a3 (£ a ]Sa - & a 2Cot)] + W 1 ( Caa) ] + Saa)2) 

3 ♦ 

- ['(w-| + h z i) (w-jCa + W 2 Sa ) + ( h x ] + *l) w 3 + V]] 

- O 2 <^(t cl sa - t c2 ca) 2 + * C 3 + 2 - 3[(* al Sa - Z a2 Ca) 2 + (h xl 


« “2 


+ *l) *i> +f EI zl V V 


3X-j 


[EI zl vf (1 -f v^)] 


= A-j - 0 


EI Z1 v i' 0 - 7 V 1 2) 


X-j = £] = 0 


(109c) 


Similarly, the differential equation for w«| is 

p 

P 1 { ( h xl + x i ) ["^(w^Ca + ajgSa) + 3 m ^a3^al Ca + & a 2 s ° 0 ] 

+ (hy-j + v-j ) [a) 3 (w]Sa - c*) 2 Ca) - £^3 (£ a ]Sa - JIggCa)] 
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+ (h z1 + w } ) [ « 2 + 0 J 2 + 2n 2 - 3fi 2 (£ g 2 + aj)] 

- V-j (o) 1 Ca + a3 2 sa)} - [(h xl + X^) (uj| Sa - u^cu) 

+ ^ h yl + v l)( w ] Ca + w 2 sa) + Wj> ~j- n 2 {(t^Sa - & c2 ca) 2 

+ z c3 + 2 - 3 C(*al Sa ' *a2 c “ )2 + C(h xl + *]) 2 

- (h x1 + x,) 2 ] + m ] (h xl + J^)} wj +| Ely-j wjw^N 

2 

4 - — E"EIyiW| (1 - _ )] + P z j = 0 , 0 < X*j < £*j (110a) 

3X-j ^ 

where must satisfy the boundary conditions 

w> ] (0) = 0 , w-j (0) = 0 (110b) 

m l {(h xl + [-W 3 (o>] Ca + u^Sa) + 3Q Z &3 (Ag-jCa + £ a2 Sa)] 

+ ^yl + V l^ Ctu 3 (o3i Sa - u> 2 Ca) - 3ft 2 (£ a jSa - fl^Cct)] 

+ (h zl + W-j) [4 + a| + 2n 2 - 3ft 2 (£ a 2 + £ a |)] 

- V, (id] C ct + 0 > 2 Sci) - — [(h xl + £]) (wjSa - ta 2 c a ) + (h yl + V, ) (u-jCa 
+ u) 2 So) + W}] - a Z ^()l c iSa - S. c2 Ca) 2 + ij + 2 - 3[(* al S« 

- «a2 c “) 2 + l a3y ^ h x1 + *l) w p + f EI y i w^wif 2 

*3^£ EI yi"i ' (1 'i w i 2)] 


EI yl W l 0 " J W ] 2 ) 


x-j = £-j = 0 
X| = = 0 


(noc) 
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In the above 


w-j = - n(s 02 C 03 + s0-|Ce2S0 3 ) - cQ^-j + ce-jse^ 

(02 = ^(s02S03 - 50 ^ 002002 ) + s 0361 + 00-100302 (Hl a ) 

C03 = ^C0-|C02 + S0]02 + 63 

£ al = c0 2 ce 3 " S0-]S02S03 

i 

z a2 = ~( c 0 2 s9 3 + sO'j (111b) 

*a3 = 00^02 

£ C 1 = -(S 02 C 03 + S9^C02Se 3 ) 

£ C 2 = S0 2 se 3 - S0-jC0 2 C0 3 (111c) 

^c3 = C0 1 C0 2 

The equations of motion and boundary conditions associated with the booms 
2, 3, and 4 are obtained from Eqs. (109) and (110) by replacing a by 
TT-a, iT+a, and 27r-a, respectively, and, of course, changing the sub- 
scripts of v-j , and w-j accordingly. 

Following the same procedure as that used to obtain Eqs. (109) 
and (110), the equation of motion and boundary conditions for V5 are 

p 5 {(h X 5 + X 5 ) [-w-|(w 2 se + W 3 CB) + 3fi^£ a -j(£ a 2SB + ^ a 3 c ^ )] 

+ (h y5 + v 5 )[u| + 4 + 2fi 2 - 3U 2 (tj + *J)] + (h z5 + 
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Men 


2 

- w 2 cb) - 3ft (^ a3 ss - £ a2 c$)] + ^ w 3 ce + w 2 se ^5 } 

^5 t(^x5 + x 5 ) (co^CB - w 3 s$) - (h Z 5 + Wg)(w2Sf3 + oj^cB) + Vg] 

• U cl + (<L c2 cs - *C3 S6 > 2 + 2 - 3 !>ai + (* a 2 ce ‘ 

^ p 5 £< h x5 + ^) 2 ' ( h x5 + x 5) 2 3> v 5 + | ei z5 v 5 v 5 2 } 


+ C-EI zS vg (1 - | v^ 2 )] + p uK = 0 


3X 


2 5 


75 


0 < < £g 


V 5 (0) = 0 , Vg(0) = 0 


EI z5 ^ + 4 CEr z5 v 5 ^ -f^ 2 )] 


X 5 = *5 


= 0 


EI z5 v 5 0 - f v 5 2 >|x 5 = * 5 “ 0 


'l 


and those for are 

P5 {(h x5 + x 5 ) [(w 3 sa - u^CBM^SB + w 3 ce) - 3ft 2 U a3 S0 - ^a^CS ) ( ^ a2 S@ 
+ ^ a 3 c &)] - (h y5 + v 5 )[wj (u^ce - oj 3 sb) - 3ft 2 * al U a2 CB - £ a3 SB) 

+ (^z5 + w 5^ ^ W 1 + ( W 2 S ^ + U 3 C 3)^ + 2ft 2 - 3ft^ [^ a ] + (£ a2 s 3 

+ *a3 c *) 2 ] > - Vg (u 2 s 6 + (D 3 C 6 ) > 


' It p 5 C"( h x5 + x 5^l + ^y5 + V 5^“3 C ® + “2 se ^ 

- 3 ^ V fl2 U c 1 + (*c 2 c6 - *c3 s6 > 2 + 2 * 3 t> a i + (» # 2 C# * *a3 ss ) 
^ p 5 [( h x5 + *5) " ( h x5 + x 5^ 2 ^ w 5 + | EI y5 w 5 w g 2 ^ 


]} - 


(112a) 

012b) 


(112c) 


]'} * 
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+ JZ !> EI y5 w 5 (1 - | W ^ 2 )] + P z5 - 0 . 0 < x 5 < l 5 


(113a) 


w 5 (0) = 0 , wj(0) = 0 


(113b) 


? EI y5 W 5 W 5 Z + aJT C«y5 w 5 0 '| W 5 2)] 


x 5 = *5 


= 0 


EI y5 w 5 (1 - I W 5 2 ) 


x 5 = *5 


- 0 


(113c) 


The equations for boom 6 are obtained from Eqs. (112) and (113) by an 
appropriate change in subscript, and by replacing a by it+8. 

b. Perturbation solution of the equilibrium problem . 

The first problem in attempting a solution of the equations of 
motion, Eqs. (108), (109), (112), and (113), is to identify the equilibrium 
configurations. To this end, we must let all the velocities and accelera- 
tions equal to zero in these equations. This leaves us three transcendental 
equations for the rotations ej (j = 1,2,3) and twelve nonlinear differential 
equations for the elastic displacements v.. , (i - 1,2,..., 6). 

We shall consider the solution of the nonlinear equilibrium problem 
in the form 


v i0< x i> = v i00 (x i } + v i01 ( V* w i0< x P = + w i01 (x i } ’ 


i = 1 ,2,. . . ,6 


OH) 


where the third subscripts on the right side of Eqs. (114) indicate the 
solution of the linearized problem if the subscript is zero and relatively 
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small perturbations if the subscript is one. It follows that the inertia 
matrix of the deformed body can be written as 

6 6 

[j(°)] = [j(°)] 0 + [jM], = E C* < 3 T CJ f ] 0 C*i3 + , 2 OifCJPiOi] 015) 

where is the inertia matrix as if the body was entirely rigid, in 

which 


J 0110 ° A 0> J 0220 ’ B 0’ J 0330 = C 0 


°0720 ‘ J 0210 ” J 0130 * J 0310 ~ J 0230 ' J 0320 " ° 


016) 


are the moments of inertia of the rigid hub, and 


J illO = p i (h yi +h zi )dx i +in i (h yi +h zi ) 

J i220 = Iq p i C (h x i + *i ) Z ^\. ]dxi ^ [ (h xi n,. )W zi ] 
J i330 = I' p i [( h xi + x 1 ) 2 +hJ i ]dx 1 .+n, i [(h xi +t 1 ) 2 +h2 i ] 


J il20 = J i210 = -J 0 p i (h xi +x i )h yi dx i- m i( h xi + 0) h yi 

rZ j 

J il30 * J i310 = ’J 0 Pi< h xi +x i) h zi dx r m i (h xi n i)h z i 

J i230 ’ J i320 = p i h yi h Z i dx i- m i h yi h zi 


i = 1,2,... ,6 (117) 


are the moments of inertia of the appendages when in undeformed state, 
expressed in terms of local coordinates. Moreover, is the change 
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in the inertia matrix due to first order elastic displacements, which has 
the elements 


J illl " L p i^ 2h yi v iOO +V iOO +2h 2i w iOO +w iOO^ dx i 


+ m. (2 h y i v. 00 +2h z 1 w . qq+w^ oq ) 


x i " M 


J i221 = / 0 p i^ h zi w i00 +w i00^ dx i +m i^ 2h zi w i00 4w 100^ 
J i 331 = L p i ( 2h yi v iOO +v iOO^ dx i 4lT1 1 ^ 2h y i v i 00 +v i 00 ^ 


x i = H 


Xi = ii 


f^i 

J il21 “ J i211 = "J 0 p i^ h xi +x i^ v i00 dx i " m i ( h xi +x i ^ v i 
r*i 

J il31 J i311 = -J 0 p i (h xi +x l) w i00 dx i ' m 1 (h xi + V w i 

f&j 

J i231 = J i321 = -] n p i (h yi W iOO +h 2 i v ioO +V iOO w iOO )dx T 


00 


00 


Xj - 


x i ” A i 


” m i ^ h yi w i00 +h zi v i00 +v i00 w i00^ 


x i = M 


i = 1 ,2,. . . ,6 


(118) 


To linearize the algebraic equations for the angles e, (j = 1 } 2 S 3) 
we would have to assume that the angles are small. This, however, is not 
always true for an arbitrary satellite, so that 1 inearization cannot be 
justified. Fortunately, it is not difficult to solve the nonlinear 
algebraic equations for the angles 9j (j = 1,2,3) by means of Newton- 
Raphson method for the moments of inertia given. As a first iteration. 
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we insert the moments of inertia of the satellite regarded as rigid into 
the tree transcendental equations for e. (j = 7,2,3), and obtain some 

J 

preliminary values for these angles. Hence, letting all terms in Eqs. 

(108) involving time derivatives equal to zero, we obtain 

4 ^12O c0 l s0 2 c0 2 + '^220 ~ ^330^ se l ce l ^ 8 1 ” 3s 9 2 ) ^^13O 00 l s®2^®2 

f 

+ ^230^ c2e l " s ^9])(3s^ 02 - c^e-j ) = 0 (119a) 

4Jno se 2 c0 2 + 4J 1 2O S0 l ^ c2 °2 " s2q 2) " 4J 220 s2q 1 s0 2 c0 2 “ 4J *3O c0 l ( c2q 2 

0 it *^t o 

- s^0 2 ) + 8J 230 se l c9 l s9 2 cs 2 - 4J 33o c 0 i s0 2 c0 2 = 0 (119b) 

~^12O^ s20 2 ” 3 c^6 2 ) + 4(Jt| 0 “ 3 22O^ 80 l 80 2^ 8 2 + 1 2 0^ ^ ^ 1 ~ 3s 2 02) 

+ 4J 230 ce l S9 2 c9 2 " ^^so 30 ! ce l (° 2e 2 " 3s2e 2) = 0 (119c) 

where [J*]q = [s]^ [J^]q Mg in which [e] 3 is given by the first of 
Eqs. (104). Regarding the angles 6j (j = 1,2,3) as known constants, we 
can linearize Eqs. (109), (110), (112), and (113) with respect to v- # w i 
and their derivatives, and solve for the perturbed elastic displacements. 
Hence, inserting Eqs. (114) into (109), we obtain the equations for v^ 
in the form 

2 

P 1 {< h x7 + X 1 ) C( w io Ca + u 20 Sa ^“10 Sa ~ <^20 Ca ^ ” ^al0 Ca 

+ ^ a 20 Sa ) (^al 0 Sa ~ A a20 Ca ^ + ( h yi + V 100^ [( w io Ca + w 20 Scx ^ 2 + ^30 

+ ta 2 - 3Q 2 ^(£ al0 ca + ^20 Sa ) 2 + *a30/B + (h zl + WlQO^O^lO 501 
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w2qCcc) - 3ft 2 £ a30^ £ al0 Sct " £ a20 Cct )^ " v 100 p l ^xl 
+ X]) S? { ( £ c io Sa “ A c20 Ca ) + ^c30 + 2 “ 3 ^ £ alO Sa “ £ a20 Ca ) 

+ £ a30^ + v 100 ^ ^ p l I-^xl + £ 1^ “ < h xl + x "|) 2 ] 

+ m 1 (^ xl + <( £ c 10 Sa - A c20 Ca)2 + *c30 + 2 ' 3 ^ £ alO Sa 

- ^ a 2o Ccc ) 2 + £ a30-^ “ EI v 100 (120a) 


where v-jqq is subject to the boundary conditions 

W°> = 0 * V 100 (0) = 0 (120b) 

2 

m l ^*Vl + £ ] ) [( w io Ca + w 20 Sa ^ w 10 Sa ~ w 20 Cot ^ " 3fi ^a!0 Ca 

+ ^ a 20 Sa )( £ alO Sa " £ a20 Ca ^ + ^yl + v 100^ w l0 Ca + w 20 So ^ 2 
+ U)| 0 + 2 ft 2 - 3n 2 ((l a10 Ca + Jt a20 sc.) 2 + Jljo) ] + (h z1 

2 

+ w T qo ) [w 3 q ( w] 0 s« - ^gCat) - 3ft £ a 30( £ al0 Sa “ ^ a 20 Cot ^ 


- v 100 fiZ ( h xl + »1> (^clO Sa - *c20 Ca > 2 + *c30 + 2 - 3 [<*alO s “ 

- ^ a 20 Ca ) 2+ l a30y J + EI v 


•M 

100 


EI v 


100 


X 1 “ 


X 1 = £ 1 


- 0 


= o 


(120c) 


The quantities wjq, £ a j 0 and £ c j 0 (j = 1 ,2*3) appearing in Eqs. (120) are to 
be calculated by using 9j (j = 1,2,3) as given by Eqs. (119). Note that 
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now primes designate total derivatives with respect to because v^qq 
depends on x alone and not on t. Similarly, for w^qq we have 

P 1 ^ h xl + x ] )C~ ul 3o( 0) lo Ca + w 20 sa ) + 3 ^ 2 *a3o( £ alO CQl + £ a20 Sa ^ 

+ (h yl + V 100 ) [u 30 (u 10 Sa - w 2 oCa) - 3n 2 a a30 (a a i 0 Sa - H a 20 Ca )] 

+ ( h z1 + w l 00 ^ K“10 + “20 ) + 2fl2 " 3fi2 ( £ al0 + 4 a20^^ 

- f>i w 100 (h x1 + x^ fl 2 {(^ c io Sot ‘ £ c20 c “) 2 + £ c 30 + 2 " 3 C( £ a10 s “ 

’ £ a20 c “) 2 + £ a30^ } + w 100 ^ { \ p l K h xl + £ 1> 2 " < h xl + x p 2] 

+ m 1 (h xl + *,)} {U clo sa - n C 2o c “) 2 + £ c 30 + 2 " 3 C( £ alO Sa 

- £ a 20 c “) 2 + £ a30^ " EI y1 w 100 = 0 (121a) 


where w-|qq is subject to the boundary conditions 

w 100 (0) = 0 , Wl ' 00 (0) = 0 0 21b) 

o 

ra l « h xl + £ P t ”“30^70°“ + w 20 sa) + 3sr£ Q 3U a ioCa + ^ a 20 Sa )] 


+ ( h y7 + v 100^ £“3oHo Sa “ Q3 20 cot ) “ 3fi2 ^a3( £ alO Sa ” ^a20 Ca )] 
+ (b zl + w 10Q ) [w 2 Q + o| 0 + 2U 2 - 3ft 2 U a f 0 + «- a 20)l 
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The differential equations and boundary conditions for v.qq and w^qq 
( i = 2,3,4) are obtained from Eqs. (120) and (121) by replacing the sub- 
scripts of and by the appropriate ones and the angle a by n-a, 
7r+a, and 27T-a, respectively. On the other hand, the differential 
equations and boundary conditions for v 5 qq and w 5 q 0 are obtained from 
Eqs. (112) and (113) in the form 

P 5 { ( h x5 + x 5^ [“^o^20 sp + “30 ce ) + ^ a io^a20 sfi + *a30 ce ^ 

+ ( h y5 + v 500^ ^20 + “30 + “ 3q2 ^a20 + ^a30^ 

+ (h z5 + w 500 ) [to 10 (a) 3 0 s$ - O) 20 ce) - 3ft 2 £ a]0 (£ a30 se - *a20 ce ^ } 


(122a) 

(122b) 

(122c) 


' V500 P5 ^ (b x5 + x 5 ) U c f 0 + U c 20 c3 - ! - C 3Q se ''> 2 
+ 2 - 3[z a 2 Q + U a20 c S ' ^a30 SB ^ 2 ^ } + v 500 ^ p 5 ^ h x5 
+ l 5 ) Z - (h x5 + x 5)^ } { * c 10 + (£ c20 ce ' *c30 se ) 2 + 2 

- 3U a10 + (£ a2 o c S - * a 30 s ®^ ' EI z 5 v 500 " 0 

where v 500 satisfies the boundary conditions 


v 50q(^ ~ ® “ 0 


500 


EI z5 v 500 


X 5 = *5 


0 ’ EI z 5 v 500 


= 0 


x 5 ^5 


as well as 

p 5 Uh x5 + x 5 ) t(^3o 5 ^ “ W20 c6 ^“20 se + - w 3 o c 3) - 3fi2 ( £ a30 s& 
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£ a20 c ^ £ a20 se + £ a30 c ^ " ^yg + v 500^ ^10 ^20^ " w 30 s ^ 

- 3fi 2 A al0 (£ a20 ce - £ a3 gS3) + (h Z 5 + w 500 ) {oj 2 q + (aj 2 gS$ + w 3Q c3) 2 

+ 2ft 2 - 3ft 2 uj 0 + (A a20 s$ + A a30 ce) 2 ]} - w 5 00 p 5 fi 2 (h x5 

+ x 5^ ^ £ clO + ( £ c20 ce ” £ c 30 s& ^ + 2 3 t £ al0 + ( £ a20 c6 ~ £ a30 S ^ ^ 

+ w^ 00 ft 2 {ip 5 [(h x5 + * 5 ) 2 - (h x5 + x 5 ) 2 ]} {£ c f 0 + U c2 o c * 

" £ c30 se ^ 2 + 2 " 3 C £ af 0 + ( £ a20 cB “ £ a30 s6 ^ 2 ^ } " EI z5 w 500 = 0 ( 123a ) 

W 0 * = 0 * W 500 (0) ““ ° (123b) 


El 


w 


y5 "500 


= 0 


El 


X 5 " *5 


w 

y5 500 


= 0 


(123c) 


x 5 = "5 


The differential equations and boundary conditions for Vggg and w^qq are 
obtained by replacing in Eqs. (122) and (123) v 5 qq and w 50Q by v 60 q and 
w 600 anc * ^ by 1T+ ^» respectively. 

On the other hand, the boundary-value problem for the perturbation 
v-jq^ is defined by the differential equation 

2 2 2 2 2 
p l v 101 fHo 0 * + w 20 Sa ^ + w 30 + ^ ~ [( £ alO Ca + £ a20 Sa ) + £ a30^ 5 

O 

+ Pl w 101 ^30^1 0 s 01 “ “20 Ca ) “ £ a30 (£ a!0 Sa ~ £ a20 Ca ^ 

+ vj Q1 {-p^ 2 (h xl + X|) ^U c -jo Sa - £ c20 Ca ) + £ c30 + 2 
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and that for the perturbation by 

P 1 V 101 ^“30 ( w lo So1 " w 20 Cct ^ ’ ^aSO^alO 501 " 

+ p l w i 0 i [“10 + “20 + 2q2 - ZqZ [ £ afo + £ a20^ 

+ wj 01 {- P] a z (h xl + x,) (U cl0 sa - n c20 c o ) + dj 0 + 2 

- 3 [U a lO s “ - £ a 20 c “) 2 + K a30^ + EI yl [ 10 w 100 w 100 

+ 5 w foo w ioo ]} + w ioi {!j2 C pi [< h xi + *i) 2 - 0\i + x i) 2 J 

+ (h x i + &i ^(^clO 501 “ ^c20 c °) + ^c30 + 2 ^[(*'a'IO Sa 


- 

s, a20 Ca ) 

+ ^30^) + 

EI yl [ 

4 Moo) 2 

+ 10 w] 00 w^ 00 

]} 


+ 

Ely^lOl 

{10 w-j 00 w^ 0C 

)) 





- 

EI yl W 101 

[1 - f (w-J 00 ) 2 D = 

- Ely] [f 

MW 3 + 10 

w i oo w i oo w i"oo 


+ 

1 W 1 00^1 00^ 





(125a) 



w 101 <°) 

= 0 

* W°> 

= 0 


(125b) 

101 

^30^1 0 s a "* ^20^°^ " 

3 ^ 2 *a30^al0 Sa " 

t a2 o c a)] 



+ 

Vioi C 

“10 + “20 + 

2q 2 - 

3^ (^alO 

+ 4 alo^ 


I 

i 

+ 

w-j Q-j {rnft 

< h xl + *l) 

( [ £ clO s “ " *c20 Ca) + “Jo + 2 

3[(£ a10 Sa 


- 

£ a 2o Ca ) 

+ ; 'a30-i ) - 

5 £J 

2 yi 

[ 2w ioo w i"oo 

+ f w ioo^=l 
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"“‘“wo , aoe 


Bu * sot niMBD 



~ EI yl w 101 ( 5w ioO w 10(P 


+ EI yl W 101 t 1 ‘ f = I EI ^ [( w inn) 2 vT 


2 yi LV "ioo' "100 


+ 2w 1 , 00 (wy 00 ) 2 ] 


x 5 = *5 


EI yi ^ ' f (-io) 2 3 wf 01 


025c) 


^ W 100 W 100 W 101^ “ ■=■ EI V/1 ( w !nn) w i ( 


2 LX yl vw 100' "100 


x 5 55 ^5 


with companion equations for and w 101 (i = 2,3,4). In a like manner 

P5V501 [ “20 + “30 + ZaZ - 3£j2 <*a20 + 2 a30 )] 

2 

+ P 5 W 501 ^10^30 s ^ “ w 20 c ^ " ^a!0^a30 s ^ " ^a20 c ^ 

+ v 501 { ‘ p 5 n2 ( h x5 + x 5> ( 2 C 10 + ( £ c20 cb ' ^o 5 * 3 ) 2 + 2 

■^alO + ( 2 a20 CB " 2 a30 S6 ^ 0 + EI z5 I- 10 v 500 v 500 + 3 v 500 v 500^ 

+ V 501 {fl2 (? p 5 ^ h x5 + 2 5) 2 ‘ ( h x5 + x 5) 2 ^) ( 2 C 10 + ( 2 c20 c6 
- ^ C30 S6) 2 + 2 - 3[£ a 2 Q + (£ a20 C6 - £ a30 sB) 2 ]) + EI 25 [If (vg 0Q ) 2 

+ 10 v 500 v 50oI } + EI z 5 v 501 ( 10 v 500 v 500^ 

' EI z5 v 501 0 - f < V 500) 2 1 ■ - ^5 Cf (^ Q0 ) 3 

+ 10 v 500 v 500 v 500 + f v 500 (v 50o) 2 ] (126a) 

V 501 ^ 8 V 5Q].^^ = ^ (126b) 
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EI z 5 f v 500 v 500 + ( v 500^ v 501 " 5 v 500 v 500 v 501 

+ n -|(V 5 00 ) 2 ] v^} 

= f EI z5 ^ v 500^ 2 v 500 + 2 v 500 ( v 500^ 


- EI z5 {[(1 - | ( V 50 O ) 2 ] v 501 


" 5 v 500 v 500 v 501 } " " " EI -c ( y inn) 


2 w ii 


2 tA z5 ' v 500' 500 J 


f at x 5 = *5 


P 5 V 501 C"lo(“20 c6 ' “’30 s8 ) ' 3s ^ Jl a10( J1 a20 c8 ‘ %30 s8 )] 

+ P 5 W 501 (“10 + (“JO 38 + “30 c8 ^ + 2 ^ 2 " ^ 2 t^alO + (^a20 S8 
+ «-a30 c8 > 2 }) + w 501 { " p 5 q2 ( h x5 + x 5> (* c 10 + (*c20 CfJ 
- *c30 s6 ) 2 + 2 ‘ 3 ^afo + (*a20 c8 - *a30SB) 2 ]/ 

+ Elyg CIO wg'oQWgQo + 5 W5Q0W500]} 

+ Wg 0 i {si 2 <^p 5 [(h x5 + i 5 ) z - (h x5 + x 5 ) 2 ]) ( a c 2 0 + U c20 ce 

" £ c30 s ^ 2 + 2 ^alO + (^a20 c8 ' ^aSO 38 ^ ^ ^ EI y5 ^ 2 ^ w 500 
+ 10 w 500 w 500^ + EI y5 w 501 ( 10 w 500 w 500^ 

' EI y5 w 501 t 1 * | Koo^ " - EI y5 tf Koo^ 

,* i 11 5 11 11 / 1 

+ 10 w 500 w 500 w 500 + — w^qq (w 500 ) ] 

w 50l(°) = 0 » w 501<°) = 0 


(126c) 


2 

i 

(127a) 

(127b) 
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EI y5 { " I ^ w 500 w 500 + f w 500^ 2 -^ w 501 " 5 w 500 w 500 w 501 


+ t 1 - | ( w 500 )2] w 501 } 


= I EI y5 [( w 5oq) 2 w 500 + 2 w 500 ( w 500^ 2 ^ 


- EI y5 {[(1 - | (w£ 00 ) 2 ] w^ Q1 


" 5 w 500 w 500 w 501 > ” ‘ f EI y5 ( w W 2 w 500 J 


| at x- = £5 (127c) 


It should be pointed out that this particular perturbation scheme 
enables us to solve first Eqs. (119) for the first approximation rotations 
9j00 (j = 1*2,3) independently of the elastic displacements. The rotations 
are then introduced into Eqs. (120) - (123), yielding the first approxi- 
mation for the elastic displacements v i0 Q and w ioo ^ = 1*2*..., 6) 
independently of one another. Inserting the first approximation v^qq 
and w-jQQ into Eqs. (124) - (127), we can obtain the corrections v^-j and 
w.jQi to the elastic displacements. The sums of these solutions yield v^g 
and w^q (i = 1,2,..., 6) according to Eqs. (114). Then, inserting v^g and 
w i0 (i = 1*2,..., 6) back into Eqs. (119), we obtain the angles e jq 
( j = 1,2,3). In the vast majority of cases, this approximation is suf- 
ficient. If not, having the new angles, we can iterate once more to 
improve the elastic displacements v iQ and w as well as the angles ejg. 


c. Liapunov stability analysis and the eigenvalue problem . 

The values 6jg, v^g, and w^g obtained above, together with sets 
of admissible functions £j(x-,0 and ^j(Xj), are subsequently introduced 
into Eqs. (49), (51), and (53), to obtain the coefficients m-^, f.^, and 
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kj k . The coefficients nij k and kj k yield the symmetric matrices [m] and 
[k], whereas using Eq. (57) the coefficients fj k yield the skew symmetric 
matrix [g]. 

From Sec. 6, if [k] represents a positive definite matrix, then 
the nontrivial equilibrium is asympototically stable. On the other hand, 
to obtain the natural frequencies, we must solve the eigenvalue problem 
in the form (66). However, before the nontrivial equilibrium can be 
determined, the system stability tested, and the natural frequencies cal- 
culated, it is desirable to use specific values for the system parameters. 
This is done in the next section. 


d. The shortening of the projections effect 

As indicated in Sec. 2, the booms are assumed to be inextensional , 
so that there is no longitudinal vibration. However, because of the 
transverse displacements, there is a shortening of the projection on the 
nominal axis of any element of length of the boom. In fact, from 
Eq. (14), the change in length of projection of any element of length dx^ 
is 


du . = - 1 
1 2 


U 


9v i 


axi 


2 r 


9w i 


3X. 


21 


dx.. , i = 1 ,2,. . . ,6 


(128) 


We shall treat this shortening as a perturbation of the spatial 
coordinate x^ , so that we can write 


X. = x i0 + X il . 


0 £ X^ q ^_&-j 9 1 = 1,2,...,! 


(129) 


where x i0 are the original spatial coordinates and x.-j are the 


70 



perturbations. From Eqs. (128), however, we conclude that the shortening 
is a second-order effect. Hence, it will not affect Eqs. (120) - (123) 
except that are to be regarded in these equations as x^g (i = 1,2,. ..,6). 
This enables us to solve for v-gg and w.gg and write the shortening of 
the projections in the form 


PiO 


"il 


du.j = 


r x iO 


[ 3v i00 

2 

-L 

3w i00 

2-r 

.0 


“il 

T 




d Si 


i = 1,2,. ..,6 


(130) 


where is a dummy variable. On the other hand, the perturbation equa- 
tions, Eqs. (124) - (127), must be modified to account for the shortening 
effect. For example, the boundary-value problem for v^g-j becomes 

p l v 101 ^lO 001 + “20 Sa > + w 30 + ^ ~ ^ 10^-j 0 ?a * ^*20^ + 

+ Pl w ioi t w 30^ w 10 Sa " w 20 Cct ^ r 3 ^a30 ^a!0 Sa ~ A a20 Ca ^ 

+ v| 01 {- P^a 2 (h xl + x 10 ) ((^10 S “ r l cZ 0 co > + *C§0 + 2 

- 3[(«. al0 sa - «. a20 Cci) 2 + i a 30 ^+ EI zl (10 v^ogvioo + 5 v^ 0 v^ 00 )} 

+ v l'oi {fi2 (|- p l C^xl + V 2 ' (h xl + x 10 )2] + m l (h xl 
+ * 10 )) (U clo s« - *c20 c “^ + ^clo + 2 ‘ 3 t(*'alO S01 ‘ £ a20 c °) 2 + s -a3o3) > 
+ EI zl C 1 ! (vf 00 ) 2 + 10 v^v^]} + EI zl 00 vi 00 vf 00 ) 1 

- Vfij, [1 - | (vi 00 ) 2 ]> = -EI zl [| (vi' 00 ) 2 + 10 v-jggVlfoovj^o 
+ | v 100 ^ v 100^ 2 ^ + J P 1 C(“io Cct + u 20 s °)(“10 Sa ‘ u 20 Ca ^ 


71 



? r 

3^ ^ a ]o Ca + ^a20 Sa ) ^alO Sa “ *a20 Ca ) j 


“10 


3v 10 0 
H- 


} 2 


faw 

dW 100 

r 

l 8? i J 



d£-j 3 0 < x-|q £ £ 1q (131) 


subject to boundary conditions (125b) and (125c), where in the latter 
must be replaced by the shortened length 

12. The RAE/B Satellite. Numerical Results . 

The general formulation of Sec. II has been used to obtain the 
nontrivial equilibrium configuration of the RAE/B satellite, to test the 
stability of equilibrium, and to calculate the natural frequencies of 
oscillation about the nontrivial equilibrium. The system parameters are 
as follows: 

Aq = 87.74 slug ft 2 , B' Q = 83.74 slug ft 2 , C Q = 18 slug ft 2 
p| = p2 = P3 = P4 - 4.348 * 10“ 4 slug ft”^ P5 ~ P6 = 4.596 
x 10" 4 slug ft“1 


m i ■ 

m 2 = m 3 = 1 TI 4 

= 2.40 * 

10” 2 slug. 

m 5 = m 6 = 

0 


4 i = 

l Z = *3 = H 

= 600 ft. 

= £ 6 = 

315 ft 



EI yi 

- EI zl 

EI y 2 

= . . . = 

EI z4 « 15. 

278 lb ft 2 

, a = 

30° 

EI y5 

= EI z5 

= EI y 6 

" EI Z 6 = 

13.889 lb 

ft 2 , 6 = 

25° 


h xl = 

h x4 = 

0.973 

ft, h x 2 = 

h X 3 = 0.878 ft, h X 5 

= h x 6 

= 0 

h yl = 

- h y4 5 

= 0.705 

ft s hy 2 

= _h y3 = " 

0.760 ft, 

h y5 = 

h y 6 

h zl = 

h z 2 = 

h z3 = 1 

h z4 = h z5 

= h z6 = 0 




fi = 

4.653 

x TO -4 

rad sec" 

1 





We shall present the results of the analyses in the order listed above. 
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a. Nontrivial equilibrium 

Inserting the above data into Eqs. (116) and (117), as indicated 
in Eq. (115), and solving Eqs. (119), we obtain 

0100 “ 0*13537 rad = 7.756 deg. 

e 2 oo 3 -5.63789 x 10 -8 rad = -3.2302 x 10’ 6 deg. 

0300 = 1-37374 x 10' 6 rad = 7.87096 x 10' 5 deg. 

We note that e^Q is caused largely by the damper booms. The fact that 
the rods are not attached at the satellite mass center turns out to have 
an insignificant effect on 6joo (j - 1,2,3). 

To evaluate the elastic displacements v^oo( x -|) and w ioo^ x i^ 

(i = 1,2,... ,6). We assume the solution of Eqs. (120) - (123) in the form 


v ioo (x f> = VioW 


w i00^ x d = z . b rioV^ x i ^ 

r= I 


i = 1,2,... ,6 (132) 


where 


<f> r (x 1 *) = A r [(cos + cosh 6 r t-j)(sin e r x.} - sinh B r x-j) 


- (sin B r &j + sinh B r £-j)(cos e r x-j - cosh 8 r Xi)] (133) 


are eigenfunctions corresponding to a bar in bending with the end x^ = 0 
fixed and having a mass m^ attached at the end x^ = . The eigenvalues 
B r &i are solutions of the characteristic equation 
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0 34) 


iru 

(7 + cos 3 r £.j cosh S r £.j) = — r- (sin & r z^ cosh g r £-j 

p i Jt i 

- cos & r Zj sinh e r J£.j) 

Moreover, the amplitudes A r are such that the eigenfunctions 4> r (x-j) are 
orthonormal, i.e., they satisfy the relation 


rz. 

J Pi* r ( x i)* S ( x i) dx i + m i^r( £ i^s^i) * 6 rs (I 35 ) 

where <$ rs is the Kronecher delta. Limiting the series in (732) to two 
terms, p = 2, the first two roots of Eq. (134) and the amplitudes A r 
corresponding to i - 1,2, 3, 4 are 

e,£, = 1.85813 A, = 0.47696 slug"’ 7 ' 2 

1 1 1 1/? 036) 

B 2 i 1 = 4.65310 A 2 - 0.03789 slug ' 

In addition, the coefficients a^g, b^.Q (i = 1,2, 3, 4) are 

Table I. 


i 

a li0 

a 2i0 

b l 10 

b 2i0 

1 

-0.13656 x 10 2 

-0.98055 x 10 -1 

0.55184 

0.46385 x 10‘ 2 

2 

0.13652 x 10 2 

0.97981 x 10' 1 

0.55188 

0.46385 x 10‘ 2 

3 

-0.13652 x 10 2 

-0.97980 x 10' 1 

-0.55187 

-0.46384 x 10' 2 

4 

0.13656 x 10 2 

0.98054 x 10' 1 

-0.55184 

-0.46385 x 10" 2 


The first two roots of Eq. (134) with m^ = 0, and the amplitudes A-j and 
A 2 corresponding to i = 5 9 6 are 


74 



$l£.j ~ 1.87571 
= 4.69414 

whereas the coefficients a w ,- 

r i 


A-j = 0.63510 

slug' 1 / 2 

A 2 * 0.04899 

slug’ 1 / 2 

and b ^ are 
n 



(137) 


Table II 


i 

a liO 

a 2i0 

b liO 

b 2i0 

5 

-0.93855 x 10’ 2 

-0.12900 x 10' 3 

0.17756 

0.70688 x 10" 3 

6 

-0.93840 x 10' 2 

-0.12900 x 10' 3 

0.17756 

0.70688 x 10 -3 


It will prove of interest to list the elastic displacements of the end 
points, as calculated by means of the linearized equations. These dis- 
placements are 


v io 0 U]) = “52.205 ft, 

v 200^2^ = 52 .192 ft, 

v 300^3^ = _52 *1 91 ft > 

v 400 ( ^4 ) * 52.204 ft, 

V 500^5^ = ~ 4.8655 x 10" 2 

v 600 U 6 ) = -4-8648 x 10' 2 


w 100 Ul) = 2.1071 ft, 

w 200^ £ 2) = 2.1072 ft, 

w 300 ^ 3 ^ = " 2* 1072 ft > 
w 400 (M) = -2.1071 ft, 

ft, w 500^5^ = 0 . 92961 ft, 

ft, w 600^5^ = °* 9 2960 ft 


The above values of v * * * * -j 00^ x i ^ anc * w i00^ x i^ ^ * 1,2,..., 6) enable 
us to solve Eqs . (119) for the angles e^Q (j * 1,2,3) and Eq. (131) and 
the companion ones for the perturbations v^fx. ), w^q-j (x^ ) (i = 1,2,..., 6). 
The resulting angles are 
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9 10 = 0.19695 rad = 11.2846 deg. 

0 2o = -6.54250 x 10" 8 rad = -3.74858 x 10' 6 deg. 
e 30 = 1.37608 x 10‘ 6 rad = 7.88438 x 10' 5 deg. 

Instead of listing the perturbations v^q-j and w-jQ-j, we shall list the 
complete solutions v-jq and in the form of the series 

2 

^ a rl - ^ r ( x f ) 

i - 1,2,-. .,6 (138) 

2 

w iO< x i) = * b riV( x i) 
r=l 

where (x^. ) are still given by Eqs. (133), in which the eigenvalues 
B r £j and amplitudes A r (r = 1,2) are given by (136) and (137). The final 
results are tabulated as follows 


Table III 


i 

a li 



a 2i 


b li 

b 2i 


1 

-0.13803 

X 

10 2 

-0.86379 x 

10' 1 

0.54865 

0.46378 x 

10' 2 

2 

0.13800 x 

10 2 

0.86313 x 

10" 1 

0.54869 

0.46378 x 

10' 2 

3 

-0.13800 

X 

10 2 

-0.86313 X 

lO' 1 

-0.54869 

-0.46378 x 

10' 2 

4 

0.13803 

X 

10 2 

0.86379 x 

10- 1 

-0.54865 

-0.46378 x 

10' 2 

5 

-0.93855 

X 

10' 2 

-0.12900 x 

10" 3 

0.17756 

0.70670 x 

10‘ 3 

6 

-0.93840 

X 

10‘ 2 

-0.12900 x 

10‘ 3 

0.17756 

0.70670 x 

10' 3 
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Moreover, the final end displacements are 


v 10 (£-|) = -52.816 ft. 


w 10 (e-|) = 2.0948 ft. 

v 20 ^ £ 2 ^ = 52.803 ft. 


w 20^2^ = 2.0950 ft. 

v 30^ £ 3^ ~ -52.803 ft. 


w _q( * 3 ) = -2.0950 ft. 

v 40 ( ) = 52.816 ft. 


w, Q (£ 4 ) = -2.0948 ft, 

v 50 U 5 ) = - 4.8655 x 10 " 2 

ft. 

w 50 U 5 ) = 0.92962 ft. 

v 60^6^ = ' 4.8648 x 10 " 2 

ft. 

W 60 (£ 6 J = 0- 92962 ft > 


and we note that the nonlinear effect is virtually zero for booms 5 and 6 
The nontrivial equilibrium is depicted in Fig. 7, where only the radial 
booms are shown because the displacements of the damper booms are in- 
significant. 

b. Liapunov stability analysis 

A stability analysis using k, as given by Eq. (60), as a testing 
function has been carried out. Essentially, the analysis reduced to 
testing the matrix [k] for positive definiteness, where the elements of 
[k] are given by Eqs. (53). The numerical values of the elements for 
the particular configuration at hand are listed in the next subsection. 
The matrix was found to be positive definite, so that the equilibrium is 
asymptotically stable . 

c. Eigenvalue problem 

Using Eqs. (49) s (51), (53), and (57), in conjunction with the 
above data, we obtain the elements 
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0 


3.42090 x 10 4 
m, . -4.78845 x 10* 1 -7.67515 x 10 

1 j J 

7.67480 x 10 8.49246 x 10' 5 

0 1.12564 x 10 5 

m 2 j 3.44664 x 10 1.57032 x 10 2 

-1.57077 x 10 2 5.48304x10 

iJ -6.31727 x 10‘ 3 2.20835 x 10 4 

00 

m, . 1.77064 x 10 2 0 

^ sj 

0 -2.89735 x 10 



e IV 


-6.31727 x 10" 3 
4.78840 x 10 -1 
6.81765 x 10 

2.20835 x 10 4 
3.44664 x 10 
-3.23228 x 10 _1 

1.31824 x 10 5 
1.77064 x 10 2 
-4.94522 x 10' 1 

1.77103 x 10 2 
0 
0 

0 

0 


4.78807 x 10 -1 
7.67512 x 10 
-8.49247 x 10" 5 

3.44739 x 10 
1.57032 x 10 2 
-5.48305 x 10 

1.77103 x 10 2 
0 

2.89734 x 10 


-7.67484 x 10 
-4.78808 x 10" 1 
6.81765 x 10 

-1.57077 x 10 2 
3.44739 x 10 
3.23042 x 10" 1 

0 

1.77130 x 10 2 
4.94522 x 10’ 1 



0 


0 


0 



•4.78845 x 10' 


3.44664 x 10 


-7.67515x10 1.57032 x 10 2 


4.78840 x 10" 1 3.44664 x 10 


7.67512 x 10 


1.57032 x 10‘ 


-4.78808 x 10' 


3.44739 x 10 


0 


0 


1 .77064 x 10 2 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1.77064 x 10 2 

0 

0 

1 

0 

0 

0 

0 

0 


0 o 0 
0 1 0 
0 o 0 


1.77103 x 10 2 0 0 

0 0 1 

0 0 0 




0 


0 


0 


0 0 0 

000 

-2.89735 x 10 0 0 

0 0 0 

0 0 0 

-4.94522 x 10' 1 0 0 

0 0 0 

1 0 0 


2.89734 x 10 0 0 

0 0 0 

0 1 0 


4.94522 x 10" 1 0 0 

0 0 0 


0 


0 


1 



0 


Table V 


1.29341 

6.74958 x 10' 3 
■6.36069 x TO' 8 

■4.23315 x 10' 7 
1.14634 x 10' 2 
■4.10887 x 10’ 2 

0 

•1.24053 x 10‘ 2 
2.17121 x 10‘ 2 


6.74937 x 10* 3 
3.36711 x 10' 11 
•5.09924 x 10' 2 

1.14630 x 10‘ 2 
■2.64834 x 10" 3 
0.88070 x 10' 8 

■1.24289 x 10' 2 
■1.40034 x 10' 2 
•4.97587 x 10' 4 


-3.36689 x 10'* 
6.74931 x 10‘ 3 
-6.36071 x 10‘ 8 

2.64825 x 10' 3 
-1.14630 x 10' 2 
-4.10889 x 10‘ 2 

1.40029 x 10' 2 
1.24288 x 10" 2 
2.17121 x 10' 2 


1.24289 x 10 -2 
0 
0 

-1.40029 x 10' 2 
0 
0 


0 

0 

0 

-9.12262 x 10 
0 
0 


9.12262 x 10' 5 
0 
0 

0 

0 

0 



-6.74965 x 10' 3 

1.14635 x 10' 2 

-1.24054 x 10' 2 

0 

0 

9 6,j 

0 

9.12263 x 10’ 5 

0 

0 

0 


0 

0 

0 

0 

0 


3.36715 x 10' 11 

-2.64836 x 10' 3 

-1.40035 x 10' 2 

0 

0 

9 7,j 

-9.12263 x TO' 5 

0 

0 

0 

0 


0 

0 

0 

0 

0 


-6.74958 x 10* 3 

-1.14634 x 10' 2 

1.24053 x 10~ 2 

0 

0 

9 8 9 j 

0 

0 

0 

-9.12258 x 10~ 5 

0 


0 

0 

0 

0 

0 


-3.36711 x ICf 11 

2.64834 x 10* 3 

1.40034 x 10‘ 2 

0 

0 

9 9,j 

0 

0 

9.12258 x 10* 5 

0 

0 


0 

0 

0 

0 

0 


-6.74931 x 10* 3 

1.14630 x 10' 2 

-1.24288 x 10* 2 

0 

0 

9 10,j 

0 

0 

0 

0 

0 


-9.12259 x ICf 5 

0 

0 

0 

0 




-3.36689 x 10" 11 

2.64824 x 10" 3 

1.40028 x 10~ 2 

0 

0 


9 1 1 ,3 

0 

0 

0 

0 

9.12259 x 

10' 5 


0 

0 

0 

0 

0 



-5.09923 x 10‘ 2 

6.88069 x 10" 8 

4.97588 x 10‘ 4 

0 

0 


g 12,j 

0 

0 

0 

0 

0 



0 

0 

-7.49379 x 10" 4 

0 

0 



6.36069 x 10“ 8 

4.10887 x 10" 2 

-2.17121 x 10’ 2 

0 

0 


g 13,j 

0 

0 

0 

0 

0 



0 

7.49379 x 10' 4 

0 

0 

0 



5.09924 x 10~ 2 

-6.88070 x 10" 8 

4.97587 x 10“ 4 

0 

0 


g 14,j 

0 

0 

0 

0 

0 



0 

0 

0 

0 

7.49380 x 

10‘ 4 


6.36071 x TO -8 

4.10889 x 10" 2 

-2.17121 x 10" 2 

0 

0 


9 15,j 

0 

0 

0 

0 

0 



0 

0 

0 

-7.49380 x 10" 4 

0 





0 


0 


Die VI 

-2.98546 x 10' 9 
5.43405 x 10" 6 
-4.36481 x 10‘ 6 

1.66943 x 10‘ 2 
1.72558 x 10' 5 
-2.79445 x 10' 7 

5.21464 x 10" 2 
6.59174 x 10' 5 
-3.34538 x 10‘ 7 

6.59442 x 10‘ 5 
0 
0 

-6.66325 x 10" 6 
0 


5.43385 x 10‘ 6 
1.53640 x 10‘ 5 
1.592715 x 10" 11 

1.72626 x 10' 5 
1.35817 x 10' 4 
4.74229 x 10" 5 

6.59442 x 10" 5 
6.66134 x 10‘ 6 
-1.64684 x 10‘ 5 

3.49900 x 10‘ 6 
0 
0 

-3.60230 x 10' 8 
0 


-1.53634 x 10' 5 
-5.43384 x 10' 6 
-4.36488 x 10' 6 

-1.35856 x 10‘ 4 
1.72628 x 10' 5 
2.79513 x 10' 7 

-6.66325 x 10' 6 
6.59447 x 10' 5 
3.34444 x 10‘ 7 

-3.60230 x 10' 8 
0 
0 

4.02153 x 10' 6 
0 


0 


0 


0 



-5.43408 x 10-5 

1.72556 x lO- 5 

6.59170 x 

10- 5 

0 

0 

k 6,j 

3.49896 x 10‘ 6 

3.60229 x 10" 8 

0 


0 

0 ' 


0 

0 

0 


0 

0 


-1.53641 x 10" 5 

1.35817 x 10" 4 

6.66133 x 

ID' 6 

0 

0 

k 7,j 

3.60229 x 10" 8 

.4-02153 x 10" 6 

0 


0 

0 


0 

0 

0 


0 

0 


5.43405 x 10' 6 

1.72558 x 10‘ 5 

6.59174 x 

10‘ 5 

0 

0 

k 8,j 

0 

0 

3.49897 x 

10" 6 

3.60230 x 10" 8 

0 


0 

0 

0 


0 

0 


1.53640 x 10' 5 

1.35817 x 10' 4 

6.66134 x 

10‘ 6 

0 

0 

k 9,j 

0 

0 

3.60230 x 

io “ 8 

4.02153 x 10‘ 6 

0 


0 

0 

0 


0 

0 


-5.43384 x 10' 6 

1.72628 x 10" 5 

6.59447 x 

10- 5 

0 

0 

k 10,j 

0 

0 

0 


0 

3.49900 x 10' 6 


-3.60230 x 10" 8 

0 

0 


0 

0 





1.53634 x 

10" 5 

-1.35856 x 

10' 4 

-6.66325 x 

10' 6 

0 

0 


k nj 

0 


0 


0 


0 

-3.60230 x 10' 8 



4.02153 x 

IQ" 6 

0 


0 


0 

0 



-1.59271 x 

10- 11 

-4.74228 x 

io “ 5 

1.64685 x 

10' 5 

0 

0 


k 12 ,j 

0 


0 


0 


0 

0 



0 


3.71129 x 

10* 5 

8.10039 x 

10 _T3 

0 

0 



-4.36481 x 

10' 6 

-2.79445 x 

10“ 7 

-3.34538 x 

10“ 7 

0 

0 

00 

k 13,j 

0 


0 


0 


0 

0 

cn 


0 


8.10039 x 

10' 13 

3.88946 x 

10* 6 

0 

0 



1.59271 x 

10- 11 

4.74229 x 

10 -5 

-1.64684 x 

10‘ 5 

0 

0 


k 14,j 

0 


0 


0 


0 

0 



0 


0 


0 


3.71129 x 10‘ 5 

-8.10040 x 10' 13 



-4.36488 x 

10‘ 6 

2.79513 x 

io ” 7 

3.34444 x 

10- 7 

0 

0 


k 15,j 

0 


0 


0 


0 

0 



0 


0 


0 


-8.10040 x 10" 13 

3.88945 x 10' 6 



Note that the elastic displacements were represented by one mode each. 

Using the formulation of Section 7, we obtain the following natural 
frequencies 


Table VII 


“ 1 ” -1 
rad sec 


1.514583 

X 

10' 2 

-1.514583 x 

10-2 

9.554415 

X 

10- 3 

-9.554415 

X 

10* 3 

9.900792 

X 

10" 3 

-9.900792 

X 

10' 3 

6.143208 x 

10" 2 

-6.143208 

X 

1CT 3 

4.319115 

X 

10' 3 

-4.319115 

X 

ID’ 3 

2.020401 

X 

10' 3 

-2.020401 

X 

10‘ 3 

2.020401 

X 

10' 3 

-2.020401 

X 

10' 3 

1.975175 

X 

NT 3 

-1.975175 

X 

lO" 3 

1.955743 

X 

10-3 

-1.955743 

X 

10" 3 

1.871461 

X 

10‘ 3 

-1.871461 

X 

10-3 

1.856504 

X 

10' 3 

-1.856504 

X 

lO’ 3 

1.856504 

X 

10' 3 

-1.856504 

X 

10* 3 

8.708815 

X 

10' 4 

-8.708815 

X 

10- 4 

6.155318 

X 

Kf 4 

-6.155318 x 

10" 4 

3.439481 

X 

10' 4 

-3.439481 

X 

10" 4 


d. Parametric study 

The stability analysis was carried one step farther by varying the 
angle a. It was found that the system was asymptotically stable for 
a = 50° , but became unstable for a = 51°. The results can be easily 
explained by the fact that in the absence of damper booms and for com- 
pletely rigid radial booms the system becomes unstable around a = 45°. 

The gravitational and centrifugal effects tend to deform the flexible 
booms in a manner that the moments of inertia about the local vertical and 
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about an axis tangent to the orbit are the same for an angle a such that 
50° < a < 51°. It should be mentioned that instability in both cases can 
be traced to angle e 3 , which tends to become large when the moment of 
inertia about the local vertical becomes larger than that about the axis 
tangent to the orbit, as at this point the "least moment of inertia" 
criterion is violated. 

The same parametric study was undertaken with respect to the 
natural frequencies. In terms of natural frequencies, instability occurs 
when at least one natural frequency (we recall that in our case the natural 
frequencies occur in pairs) reduces to zero. Here again the system be- 
comes unstable for 50° < a < 51°, thus corroborating the results obtained 
by the Liapunov stability analysis. 

13. Summary and Conclusions 

Two new theories for studying the motion characteristics of a 
rotating system with flexible parts about undeformed equilibrium have been 
developed. The first is qualitative and the second quantitative. 
Specifically, the first represents a stability theory and the second a 
method for obtaining the system natural frequencies. 

The stability theory is based on the Liapunov direct method and 
makes use of modal analysis to represent elastic displacements. The 
novelty of the formulation lies in the fact that for the first time a 
nontrivial equilibrium is considered in conjunction with the Liapunov 
direct method for a stability analysis of spinning flexible bodies capable 
of large deformations. 

The stability analysis can be divided into two major parts: the 

evaluation of the nontrivial equilibrium and the stability analysis itself. 
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When the body is capable of large deformations, nonlinear algebraic 
and differential equations must be solved for the rotational and elastic 
displacements, respectively, where these displacements define the 
equilibrium configurations of the system. Because the problem is one of 
stability about nontrivial equilibrium, it is necessary to expand the 
Liapunov function about that equilibrium. Assuming small displacements 
from equilibrium, the problem reduces to the evaluation of a Hessian 
matrix at the nontrivial equilibrium and testing the matrix for sign 
definiteness by means of Sylvester's criterion. It should be pointed out 
that the size of the Hessian matrix depends on the number of eigenfunctions 
used to represent the elastic displacements. 

The method for obtaining the natural frequencies of the system 
makes use of the variational equations about the nontrivial equilibrium. 
Then the set of second-order differential equations is converted into a 
set of twice the number of first-order differential equations. The 
associated eigenvalue problem yields the system natural frequencies. 

The two methods are quite general in scope, and can be used for 
testing stability and calculating the natural frequencies of a large 
variety of hybrid systems. As an application, the theory has been used to 
test the stability of the RAE/B satellite. First, the nonlinear equations 
have been solved for the nontrivial equilibrium configuration, and then 
this configuration has been used to evaluate the associated Hessian 
matrix. The satellite was found to be stable. Then one of the systems 
parameters has been varied to predict at which point the equilibrium be- 
comes unstable. The results are in line with the expectations. In 
addition, the system natural frequencies for oscillation about the 
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deformed equilibrium were calculated. The parametric study used in con- 
junction with the Liapunov stability analysis was used to examine how the 
frequencies are affected. The study resulted in the same instability 
statement. 
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FIGURE 2 - ELASTIC DISPLACEMENTS 



FIGURE 3- 0R8ITAL AXES AND BODY AXES 
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